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CHAPTER  1 
INTRODUCTION 


1.1  BACKGROUND 

The  Electromagnetic  Compatibility  Analysis  Center  (ECAC)  is  responsible  for  developing 
capabilities  for  analyzing  the  compatibility  of  communications  and  electronic  equipment  that  make 
use  of  the  electromagnetic  spectrum.  This  responsibility  includes  the  development  of  mathematical 
models  of  propagation  phenomena  for  a  wide  range  of  applications.  Some  of  the  earliest  computer- 
automated  models  developed  at  ECAC  were  purely  theoretical  techniques  for  calculating  radio-wave 
propagation  loss  between  antennas  above  an  imperfectly-conducting  spherical  earth.  Later,  a  model 
was  developed  that  required  a  point-by-point  representation  of  the  terrain  elevations  at  intervals 
along  the  irregular  profile  of  the  earth  in  the  great-circle  plane  between  the  antennas.  Both  types  of 
models,  whether  for  use  over  smooth  spherical  earth  or  over  rough  earth,  must  evaluate  the 
geometry  of  the  path,  determine  the  appropriate  mode(s)  of  propagation,  and  apply  suitable 
algorithms  to  calculate  the  path  loss  for  the  dominant  mode. 

1.1.1  Spherical-Earth  Models 

In  a  spherical-earth  model,  the  terrain  profile  is  represented  as  the  arc  of  a  circle  whose  radius 
may  depend  on  the  gradient  of  atmospheric  refractivity.  Such  a  model  is  valid  for  propagation  over 
a  smooth  sea  or  lake  or  even  over  a  relatively  smooth  plain  if  the  frequency  is  low  enough  that  the 
wavelength  is  long  compared  with  the  dimensions  of  the  terrain  irregularities.  One  of  the  earliest 
ECAC  spherical-earth  models  is  the  Smooth  Curve  Smooth  Earth  (SCSE)  model.  Later,  a 
simplified  spherical-earth  model,  intended  for  lower  frequencies  and  antenna  heights,  and  a 
modification  of  a  simple  method  for  predicting  tropospheric  scatter  propagation  loss  were  combined 
with  SCSE  in  the  Integrated  Propagation  System  (IPS).'  '  At  about  the  same  time  a  theoretical 
ground-wave  model  called  NLAMBDA,  capable  of  accommodating  antennas  many  wavelengths  in 


'-'Meidenbauer,  R.,  Chang,  S.,  and  Duncan,  M.,  A  Status  Report  on  the  Integrated  Propagation  System 
(IPSL  ECAC-TN-78-023,  DoD  ECAC,  Annapolis,  MD,  October  1978. 
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height,  was  developed.  Later,  the  NLAMBDA  capability  for  calculation  of  path  loss  was  used  as 
the  foundation  for  the  current  ECAC  standard  for  spherical-earth  models,  the  ECAC  Ear-Field 
Smooth-Earth  Coupling  Code  (EFFSECC).*"^  EFFSECC  combines  the  theoretical  path-loss 
calculation  for  an  ideal  spherical  earth  with  the  gain  and  radiation  characteristics  for  a  package  of 
commonly  used  linear  antennas. 

NLAMBDA  and  EFFSECC  provide  accurate  theoretical  predictions  of  ground-wave  loss  or 
electric  field  strength  between  antennas  on  a  spherical  earth.  However,  for  many  applications 
involving  the  analysis  of  large  environments  of  equipment  and  the  estimation  of  the  necessary 
separation  distance  required  to  prevent  interference  between  systems,  the  IPS  was  used.  The  IPS 
yields  acceptable  accuracy  in  situations  requiring  conservation  of  computation  time. 

In  1990,  the  Spherical-Earth  Model  (SEM)  was  designed  to  replace  the  IPS  for  calculation  of 
propagation  loss  over  smooth  homogeneous  terrain.  The  model  design  is  considerably  more 
compact  than  IPS  (line  count  is  reduced  by  80  percent),  yet  the  program  is  thoroughly  commented 
(annotated  with  explanatory  comments)  for  easy  integration  into  analysis  programs.  The  modular 
construction  of  the  model  makes  it  possible  to  select  subroutines  for  specific  applications.  Many  of 
the  modular  SEM  subroutines  were  specifically  designed  for  both  SEM  and  the  rough-earth  model 
described  below  in  paragraph  1.1.2. 

1.1.2  Rough-Earth  Models 

Rough-earth  models  take  into  account  the  effects  of  the  irregular  terrain  along  the  propagation 
path  and  in  the  vicinity  of  the  antennas.  Thus,  rough-earth  models  require  information  regarding 
terrain  elevations  in  the  region  of  interest.  The  development  of  rough-earth  models  began  in  the 
early  1960s  with  the  formulation  of  techniques  for  extracting  the  necessary  parameters  for  evaluation 
of  the  path  from  great-circle  terrain  profiles,  selection  of  the  appropriate  mode  of  propagation,  and 


’■^Maiuzzo,  M.  A.  and  Frazier,  W.  E.,  A  Theoretical  Ground  Wave  Propagation  Model  -  NLAMBDA 
Model.  ESD-TR-68-315,  DoD  ECAC,  Aimapolis,  MD,  December  1968. 

‘■^Marcus,  S.,  Description  of  the  ECAC  Far-Field  Smooth-Earth  Coupling  Code  (EFFSECC). 
ECAC-TN-80-014,  DoD  ECAC,  Armapolis,  MD,  August  1980. 
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calculation  of  the  loss.  Many  of  the  concepts  and  algorithms  employed  in  the  rough-earth  model 
known  as  the  Terrain-Integrated  Rough-Earth  Model  (TIREM)  were  based  on  work  done  at  the 
Central  Radio  Propagation  Laboratory  (CRPL)  of  the  National  Bureau  of  Standards  (NBS).^'^  CRPL 
has  since  become  the  Institute  for  Telecommunication  Sciences  (ITS)  of  the  National 
Telecommunications  and  Information  Administration  (NTIA). 

TIREM  provided  the  first  means  for  determining  terrain-dependent  propagation  path  loss  entirely 
by  automated  computer  processes.  Early  TIREM  versions  included  three  line-of-sight  (LOS)  modes 
(including  free  space)  and  three  beyond  LOS  (BLOS)  modes.  There  was  also  a  spherical-earth 
mode  which  served  as  a  recovery  model  when  rough-earth  engineering  models  were  not  applicable. 

Over  the  ensuing  years,  the  TIREM  modeling  was  revised  with  new  capabilities  and  features. 

By  1983,  the  model  recognized  unique  path  geometry  situations  requiring  loss  algorithms  for  three 
LOS  modes,  five  diffraction  modes,  and  five  combination  modes  involving  tropospheric  scatter  or 
diffraction.  This  TIREM  version  could  calculate  the  basic  transmission  loss  over  irregular  terrain 
for  frequencies  between  40  and  20,000  MHz.  Comparisons  with  propagation  loss  measurements 
made  by  ITS  showed  root  mean  squared  (RMS)  prediction  errors  for  various  measurement  regions 
of  10  to  43  dB.  Later  in  1983,  TIREM  was  modified  to  allow  a  low  frequency  limit  of  20  MHz  and 
a  low  antenna  height  limit  of  0.5  m. 

In  response  to  observations  that  several  of  the  prediction  modes  for  the  region  BLOS  were 
inaccurate,  an  investigation  of  the  beyond  LOS  methods  was  performed.  The  BLOS  calculations 
were  improved  by  combining  the  diffraction  calculations  into  a  single  unified  mode  based  on 
diffraction  over  multiple  knife  edges.  Also,  in  the  transition  region  between  propagation  by 
diffraction  and  propagation  by  troposcatter,  a  single  new  mode  was  added  in  place  of  several 
transition  modes.  A  comparison  of  propagation  predictions  of  the  revised  version  of  TIREM  with 
measured  propagation  data  showed  that  it  significantly  outperformed  the  existing  version  of  TIREM. 


'■^Rice,  P.  L.,  Longley,  A.  G.,  Norton,  K.  A.,  and  Barsis,  A.  P.,  Transmission  Loss  Predictions  For 
Tronosnheric  Communication  Circuits.  Volumes  1  and  2.  NBS  Technical  Note  101,  CRPL,  NBS, 
Boulder,  CO,  May  1965,  Revised  January  1967. 
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This  revised  TIREM  version  became  EC  AC  TIREM  Version  1.0  and  was  described  in  the  first 
TIREM  Handbook.’-^ 

TIREM  version  1.0  sometimes  exhibited  an  unexpected  decrease  of  loss  with  increasing  distance 
for  paths  over  smooth  water.  This  discontinuity  was  caused  by  mode  switching  between  LOS  and 
rough-earth  diffraction  modes  when  the  diffraction  occurred  on  a  smooth  water  surface.  TIREM 
Version  2.0  added  a  smooth  spherical-earth  propagation  mode  for  path  segments  over  smooth  water 
surfaces.  Losses  for  water  segments  are  combined  with  multiple  knife-edge  diffraction  losses  for 
land  segments,  weighting  each  by  the  proportionate  distances  of  water  and  land  in  the  entire  path. 
These  measures  assure  that  a  correct  loss  calculation  is  made  for  diffraction  over  the  surface  of  a 
smooth  earth  and  that  a  smooth  transition  occurs. 

In  order  to  facilitate  integration  into  analysis  programs,  TIREM  Version  2.0  was  completely 
rewritten  as  a  callable  subroutine  in  standard  ANSI  FORTRAN  77  and  became  TIREM  Version 
3.0.^  '’  The  engineering  design  of  some  of  the  algorithms  was  also  revised  to  simplify  the 
calculations  and  preserve  the  model  accuracy.  For  example,  the  smooth  spherical-earth  mode 
calculations  were  improved  and  validated  with  propagation  measurements.  This  allowed  the  lower 
frequency  limit  to  be  extended  to  1  MHz.  Version  3.0  also  simplifies  the  LOS  calculations. 

Recent  efforts  to  improve  TIREM  have  concentrated  on  simplification  and  reduction  of  the 
computation  time.  Thus,  the  model  is  more  suitable  for  application  programs  requiring  many  loss 
calculations  such  as  the  analysis  of  large  environments  of  equipment  or  the  plotting  of  propagation- 
related  performance  contours.  Although  past  efforts  have  produced  a  model  with  a  near-zero  mean 
prediction  error,  the  standard  deviation  is  11  dB,  a  figure  that  may  be  unacceptable  for  some 
applications.  Further  improvement  of  the  error  statistics  of  TIREM  has  proven  to  be  difficult  using 
elementary  yet  expedient  algorithms.  EC  AC  continues  to  investigate  other  more  sophisticated 
techniques  for  propagation  prediction  such  as  the  geometric  theory  of  diffraction  (GTD),  integral 


^■^Benoit,  G.,  Terrain-Integrated  Rough-Earth  Model  (TIREM)  Handbook,  ECAC-HDBK-86-076,  DoD 
EC  AC,  Annapolis,  MD,  September  1986. 

‘■^Sciandra,  R.  M.,  TIREM/SEM  Programmer’s  Reference  Manual.  ECAC-CR-90-039,  DoD  ECAC, 
Annapolis,  MD,  July  1990. 
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equation  methods,  and  methods  based  on  the  parabolic  equation  (PE)  approximation  to  the  solution 
of  the  Helmholtz  equation  for  electromagnetic  propagation.  Some  of  these  methods  have  the 
potential  not  only  to  account  for  the  effects  of  an  inhomogeneous  atmosphere,  but  also  to  account 
for  three-dimensional  terrain  effects  outside  the  great-circle  plane. 

1.2  PURPOSE 

The  purpose  of  this  Handbook  is  to  provide  a  description  of  the  current  versions  of  TIREM  and 
SEM,  including  the  engineering  formulations  on  which  the  models  are  based.  Both  TIREM  and 
SEM  are  modular  programs  with  many  common  subroutines  that  conform  with  ANSI  FORTRAN  77 
standards.  The  two  programs  are  designed  for  easy  integration  into  other  programs.  Both  models 
have  been,  and  continue  to  be,  incorporated  into  performance  evaluation  programs  as  replacements 
for  earlier  versions  of  rough-earth  and  smooth-earth  propagation  models. 

SEM  is  considerably  faster  than  TIREM  because  it  does  not  require  the  extraction  of  a  terrain 
profile  from  elevation  data.  However,  it  is  theoretically  valid  only  when  used  for  propagation  over 
a  smooth,  spherical  earth.  It  is  commonly  used  when  approximate  loss  predictions  are  required  or 
when  terrain  elevation  data  are  unavailable. 

TIREM  and  SEM  can  be  used  as  independent  propagation  models  for  rough-earth  or  smooth- 
earth  applications,  respectively.  However,  subroutines  from  SEM  are  extensively  used  in  TIREM. 
Therefore,  the  engineering  theory  of  TIREM  and  SEM  are  described  together  in  this  Handbook. 

The  models  are  modular  and  contain  26  conunon  subroutines.  Both  models  conform  with  ANSI 
FORTRAN  77  and  use  argument  lists  rather  than  common  blocks  whenever  possible.  The  models 
described  in  this  Handbook  are  TIREM  Version  3.09  and  SEM  Version  1.07,  and  reflect  corrections 
and  enhancements  to  the  original  models. 
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1.3  OVERVffiW  OF  MODEL  OPERATION 

Figure  1-1  is  a  top-level  block  diagram  of  TIREM  and  SEM;  detailed  flowcharts  are  presented 
in  Reference  1-6.  Both  models  provide  estimates  of  the  propagation  loss  over  an  imperfectly 
conducting  earth.  Both  models  consider  the  effects  of  free-space  spreading,  reflection,  surface-wave 
propagation,  diffraction,  tropospheric-scatter  propagation,  and  atmospheric  absorption,  as 
appropriate.  Only  TIREM  evaluates  the  effects  of  irregular  terrain  along  the  great  circle  path. 

SEM  regards  the  surface  of  the  earth  as  a  perfectly  smooth,  imperfectly  conducting  sphere.  Neither 
model  considers  ducting  phenomena,  fading,  or  ionospheric  (skywave)  propagation.  Calculations  of 
performance  improvement  that  would  result  from  angle,  space,  or  frequency  diversity  are  not 
performed.  Lastly,  neither  model  explicitly  addresses  attenuation  due  to  rain,  foliage,  or  man-made 
obstacles.  However,  if  the  terrain  elevation  profile  supplied  to  TIREM  has  been  adjusted  by  adding 
the  height  of  known  man  made  obstacles  or  trees,  their  effect  on  diffraction  will  be  included  in  the 
predicted  path  loss. 

1.3.1  TIREM  Overview 

TIREM  is  designed  to  calculate  basic  median  propagation  loss  (hereafter  referred  to  as  path  loss) 
over  irregular  terrain  for  frequencies  between  1  and  20,000  MHz.  Because  it  can  consider  the 
terrain  elevations  along  the  path,  TIREM  is  one  of  the  most  frequently  used  propagation  models  at 
EC  AC.  TIREM  uses  as  input  a  terrain  profile  described  by  a  number  of  discrete  points,  the 
position  of  which  is  specified  by  a  distance  from  the  transmitter  and  an  elevation  above  mean  sea 
level.  Also  required  is  information  on  the  transmitter  (antenna  height,  frequency,  and  antenna 
polarization),  the  receiver  (antenna  height),  atmospheric  constants  (surface  refractivity  and  humidity) 
and  ground  constants  (relative  permittivity  and  conductivity).  There  is  also  a  profile  indicator  flag 
to  save  time  when  a  profile  is  an  extension  of  the  existing  profile.  A  description  of  TIREM 
parameters  and  their  valid  ranges  is  presented  in  Appendix  A. 
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Figure  1-1.  Top-Level  Block  Diagram  of  TIREM  and  SEM 
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1.3.2  SEM  Overview 

SEM  is  a  stand-alone  spherical-earth  model  that  disregards  the  effects  of  terrain  features.  The 
basis  of  SEM  is  the  spherical  earth  propagation  module  contained  in  the  subroutine  SEARTH  and  its 
ancillary  subroutines.  SEARTH  is  a  callable  subroutine  useful  for  spherical-earth  propagation 
prediction.  SEARTH  is  also  used  by  TIREM  to  calculate  propagation  loss  over  sea  or  smooth 
earth.  In  the  stand-alone  application,  SEM  inputs  are  the  distance  between  transmitter  and  receiver 
and  information  concerning  the  transmitter  (antenna  height,  frequency,  and  antenna  polarization),  the 
receiver  (antenna  height),  atmospheric  constants  (surface  refractivity  and  humidity)  and  ground 
constants  (permittivity  and  conductivity).  SEM  calculates  and  provides  as  output  the  basic 
propagation  loss  for  antenna  terminals  located  on  a  smooth  spherical  earth. 

For  a  path  within  LOS,  the  LOS  loss  relative  to  free  space  is  the  smaller  of  the  space- wave  loss 
or  the  surface  wave  loss.  For  a  path  extending  well  into  the  diffracted  region,  the  loss  is  the 
smaller  of  the  diffraction  loss  or  the  troposcatter  loss.  Within  a  region  that  spans  the  horizon, 
where  neither  the  LOS  loss  nor  the  diffracted  wave  loss  is  entirely  valid,  the  loss  is  obtained  by  an 
empirical  splice  function. 

If  the  path  is  BLOS,  the  total  diffraction  loss  is  the  sum  of  the  free-space  loss,  the  loss  relative 
to  free  space,  and,  if  the  frequency  exceeds  10  GHz,  the  loss  due  to  atmospheric  absorption.  The 
total  path  loss  is  set  to  the  smaller  of  the  total  diffraction  loss  and  the  total  tropospheric  scatter  loss. 

1.4  CONVENTIONS 

1.4.1  Unit  of  Measure  Conventions 

Throughout  this  document,  metric  units  of  measurement  are  used.  Distances  and  heights  are  in 
meters,  unless  indicated  otherwise.  All  angles  are  in  radians.  Frequency  is  in  MHz.  Propagation 
loss  (or  path  loss)  is  the  loss  between  loss-free  isotropic  antennas  whose  polarization  is  the  same  as 
the  transmitted  wave  and  is  expressed  in  decibels. 
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1.4.2  Software  Version  Conventions 

Each  computer  output  carries  the  version  number  of  the  model.  With  each  minor  change,  the 
decimal  will  be  raised  (i.e.,  in  TIREM,  3.00,  3.01,  3.02,  etc.,  and  in  SEM,  1.00,  1.01,  1.02,  etc.) 
and  with  each  major  change  the  number  will  be  reset  to  the  next  larger  whole  number  (i.e.,  in 
TIREM,  4.00,  5.00,  etc.,  and  in  SEM,  2.00,  3.00,  etc.). 
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CHAPTER  2 

COMMON  CONCEPTS  AND  TECHNIQUES 


2.1  GENERAL 

Within  TIREM  and  SEM,  there  are  a  number  of  concepts  and  techniques  that  are  identical  in 
both  programs.  Some  components  common  to  both  programs  are  free-space  propagation,  effective 
earth  radius,  molecular  absorption,  and  tropospheric  scatter  propagation.  These  concepts  and 
techniques  are  discussed  in  this  chapter.  Because  the  spherical-earth  calculations  that  are  integral  to 
SEM  are  also  part  of  TIREM,  these  algorithms  are  also  common  to  both.  However,  the 
methodology  for  calculation  of  propagation  over  spherical  earth  is  fundamental  to  SEM  and  it  is 
discussed  in  Chapter  3. 

2.2  FREE-SPACE  PROPAGATION 

The  fundamental  concept  in  propagation  formulations  is  propagation  in  free  space,  that  is,  in  a 
region  free  of  all  objects  that  might  absorb  or  reflect  electromagnetic  energy.  For  an  ideal  isotropic 
anteima,  which  radiates  uniformly  in  all  directions,  the  emitted  power  P,  is  distributed  uniformly 
over  spheres  of  increasing  radius.  The  power  density  S  at  a  distance  d  is: 


S 


Pt 

47rd^ 


(2-1) 


With  X  the  wavelength,  the  power  P^  available  to  an  isotropic  receiving  anteima  having  an  effective 
aperture  of  X^/47r  is: 


(2-2) 
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From  Equations  2-1  and  2-2,  the  free-space  transmission  ratio  is  therefore: 


Pr 

Pt 


(2-3) 


and  the  free-space  propagation  loss  Lps  between  isotropic  antennas  is: 


^Ps  =  10  log  (^)  = 


20  log  (M) 


(2-4) 


All  formulas  in  subsequent  sections  that  calculate  the  electric  field  strength  are  normalized  to  the 
free-space  electric  field  strength,  E^,,  at  the  same  point.  In  all  cases,  however,  the  final  quantity 
required  is  the  propagation  loss  between  antennas,  rather  than  the  value  of  the  electric  field  strength 
at  a  point.  This  loss  is  expressed  as  the  sum  of  the  free-space  loss  from  Equation  2-4  and  the 
excess  loss  relative  to  free  space. 

2.3  EFFECTIVE  EARTH  RADIUS 

The  properties  of  the  atmosphere  vary  with  elevation.  The  atmospheric  property  of  concern  in 
this  subsection  is  the  variation  of  refraction  index  with  elevation.  In  a  well  mixed  atmosphere,  the 
refractivity  diminishes  with  elevation  and  this  gradient  causes  the  radio  rays  to  bend.  However,  the 
curvature  of  the  ray  path  can  be  modified  by  a  simple  geometrical  transformation  to  a  change  in  the 
curvature  of  the  earth,  so  the  propagation  of  rectilinear  rays  is  effectively  restored.  This  is  the 
effective-earth  radius  concept  used  in  tropospheric  wave  propagation  (see  Figure  2-1).  The 
transformation  is  accomplished  by  substituting  an  effective  earth  radius  a^  for  the  actual  earth  radius 
a^  in  path  geometry  calculations: 


a 


e 


(2-5) 
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LINEAR 


f 


The  distance  a,,  is  the  actual  earth  radius  (=  6.372  x  10^  meters)  and  k  is  the  effective  earth  radius 
factor.  Within  approximately  3200  m  of  the  surface  of  the  earth,  the  gradient  of  atmospheric 
refractivity  normally  decreases  at  a  rate  that  is  linear  with  height.  The  slope  tends  to  be  correlated 
with  the  value  of  refractivity  measured  at  the  surface.  At  higher  altitudes,  the  gradient  decreases 
exponentially  with  height.  Thus,  TIREM  and  SEM  employ  either  of  two  forms  for  the  effective 
earth  radius  k  factor,  one  for  antennas  below  3200  m  in  the  linear  portion  of  the  atmosphere,  and 
one  for  use  if  either  antenna  extends  into  the  exponential  atmosphere.  Thus,  k  is  defined  as: 


k  = 


kn„ 


1  - 


0.04665exp(0.005577N5) 
1  +  lO'^N. 


-1 


kexp  =  4/3 


Ah, 

1  -  ■ 

Ah,  ] 

h,  +  ^hjh. 

- 1 

+ 

.if 

hj,h2<  3200m 


max(hj,h2)>3200  m 


(2-6) 


where  the  variables  are  defined  in  the  next  two  paragraphs. 
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The  factor  k,in  is  determined  by  the  nearly  constant  gradient  of  refractivity  near  the  surface  of  the 
earth  for  heights  of  less  than  3200  m.  Nj  is  the  surface  value  of  refractivity  at  the  lower  of  the 
antenna  elevations;  hj  and  h2  are  the  structural  antenna  heights  at  those  elevations. 

The  factor  kg^p  is  derived  from  the  CRPL  Exponential  Atmosphere^  *  as  follows.  In  Reference 
2-1,  the  horizon  distance  dn  for  a  given  (antenna)  height  h  is  tabulated  for  a  set  of  values.  Also 
tabulated  is  the  incremental  height  Ah  required  to  produce  the  same  horizon  distance  that  would 
result  for  a  k  =  4/3  earth-radius  factor,  so  that: 

=  ^(4/3)a„(h-Ah)  =  V2k,,pa„h 


and  hence 


k^^p  =  4/3  (1-Ah/h) 


(2-8) 


This  is  the  effective  earth  radius  factor  for  one  antenna  height.  For  two  anteimas  hj,  i  =  1,  2,  each 
earth  radius  factor  k^^pj  is  weighted  by  the  ratio  of  its  horizon  distance  d^ ;  to  the  total  horizon 
distance  d^^  +  d^  j  as  follows: 


kexp,i  =  4/3(l-AhA) 


dH.i 


(2-9) 


so  that 


Kp  =  ke,p.i+k^,p,2  =  4/3 


1-- 


Ah, 


Ah2 


(2-10) 


^■*Bean,  B.  R.  and  Thayer,  G.  D.,  CRPL  Exponential  Reference  Atmosphere.  US  Dept,  of  Commerce, 
National  Bureau  of  Standards,  NBS  Monograph  4,  October  1959. 


2-4 


ECAC-HDBK-93-076 


The  transition  height  of  3200  m  in  Equation  2-6  is  selected  such  that  at  that  height  kiin  = 
for  all  values  of  Nj. 

The  subroutine  FAKTOR  has  stored  a  subset  of  the  tabulated  Ah  data  of  Reference  2-1  (for 
initial  angle  6  =  0,  corresponding  to  the  horizon  distance)  for  Nj  =  200.0,  252.9,  289.0,  313.0, 
344.5,  377.2,  404.8,  450.0,  and  for  incremental  heights  up  to  475,000  m.  For  intermediate  values 
of  height  and  refractivity.  Ah  is  derived  by  linear  interpolation.  The  k  factor  given  in  Equation  2-6 
is  implemented  in  subroutine  RFC2AK. 

2.4  MOLECULAR  ABSORPTION 

In  a  clear  atmosphere,  oxygen  and  water  vapor  absorb  microwaves.  Oxygen  has  a  magnetic 
interaction  with  the  incident  radiation  because  the  O2  molecule  is  normally  in  a  paramagnetic  state, 
and  water  vapor  contributes  because  of  the  electric  polarity  of  the  H2O  molecule.  There  are  certain 
frequency  regions  where  the  absorption  for  both  oxygen  and  water  vapor  is  abnormally  large 
because  of  resonance.  Figure  2-2  gives  an  example  of  the  specific  attenuation  in  dB/km  of  H2O  and 
O2  versus  frequency.  The  plots  are  for  the  average  atmospheric  conditions  at  sea  level  in  temperate 
latitudes  in  summer. 

At  higher  altitudes,  the  specific  attenuation  typically  decreases  at  a  rate  that  is  dependent  on  the 
absolute  temperature  T(H),  the  pressure  P(H),  and  on  the  density  p(H)  of  the  water  vapor  at  a 
height  H  (km  above  mean  sea  level).  Consequently,  molecular  absorption  varies  as  a  function  of 
height  and  its  calculation  is  dependent  upon  the  selection  of  an  atmospheric  molecular  model. 
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PRESSURE: 

1  ATMOSPHERE 

TEMPERATURE: 

20°C 

WATER  VAPOR: 

7.5  G/M^ 

Figure  2-2.  Specific  Attenuation  of  Oxygen  and  Water  Vapor 
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A  particular  molecular  model  used  in  the  meteorological  and  atmospheric  physics  community  is 
the  US  Standard  Atmosphere,  1962^'^  as  partially  reproduced  in  Table  2-1,  where  P^,  and  are 
tabulated  with  height  H  as  a  parameter.  This  molecular  model  is  used  to  derive  the  combined  H2O 
and  O2  specific  attenuation  (absorption  coefficient)  7  which  is  the  molecular  absorption  per  unit  path 
length  at  the  height  H.  The  coefficients  7  used  in  TIREM  and  SEM  were  generated  using  the 
Millimeter  Wave  Propagation  Prediction  Model  (MMWPROP)^-^'^'^  and  stored  in  data  arrays  in 
subroutine  GAMMAR.  Absorption  coefficients  are  stored  for  discrete  values  of  frequency  (f  =  1, 
5,  10,  15,  20  GHz),  height  (H  =  0,  2,  4,  6,  8,  10,  15,  20,  25,  30  km)  and  humidity  at  sea  level 
(p(0)  =  0,  10,  30,  50,  no  g/m^).  Intermediate  values  are  linearly  interpolated  from  these 
tabulations. 


Table  2-1.  US  Standard  Atmosphere,  1962  (from  Reference  2-2) 


Height  H  (km) 

Pressure  (mbar) 

Temperature  T,  (K) 

Water  Vapor  p,  (g/mj) 

0 

1.013E+03 

288.1 

5.9E+00 

1 

8.986E+02 

281.6 

4.2E+00 

2 

7.950E+02 

275.1 

2.9E+00 

3 

7.012E+02 

268.7 

1.8E+00 

4 

6.166E+02 

262.2 

l.lE+00 

5 

5.405E+02 

255.7 

6.4E-01 

6 

4.722E+02 

249.2 

3.8E-01 

7 

4,lllE+02 

242.7 

2.1E-01 

8 

3.565E+02 

236.2 

1.2E-01 

9 

3.080E+02 

229.7 

4.6E-02 

10 

2.650E+02 

223.2 

1.8E-02 

^■^McClatchy,  R.  A.,  et  al.  Optical  Properties  of  the  Atmosphere  (Third  Edition),  AFCRL-72-0497, 
Air  Force  Cambridge  Research  Laboratories,  August  1972. 

^■^Meidenbauer,  R.,  A  Millimeter  Wave  Propagation  Prediction  Model-MMWPROP, 
ECAC-TN-81-046,  DoD  ECAC,  Annapolis,  MD,  September  1981. 

^■^Liebe,  H.  J.,  "An  Updated  Model  for  Millimeter  Wave  Propagation  in  Moist  Air,"  Radio  Science. 
Vol.  20,  No.  5,  1985,  pp.  1069-1089. 
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The  humidity  at  sea  level  p(0)  is  obtained  from  the  local  absolute  humidity  Plocal  (in  g/m^ 
required  as  input)  from: 


P(0) 


1.7  X  10^  ^LOCAL 


Ps 


(2-11) 


where  the  temperature  and  the  vapor  density  p^  are  again  dependent  upon  the  height  H  as 
tabulated  in  Table  2-1.  T^  and  Ps  are  approximated  from  the  tabulated  values  by: 

Tg  =  -6.48  H  +  288.1  (2-12) 


and 


5.9  (0.7011)”  0<H<2  (2-13) 

7.6455  (0.6159)”  H>2 


Because  the  molecular  absorption  coefficient  varies  with  height  along  the  radio  ray  path  due  to 
variations  in  atmospheric  temperature,  pressure,  and  water  vapor  content,  the  total  atmospheric 
absorption  is  obtained  by  integration  along  the  entire  path.  Since  the  height  is,  for  a  particular  path, 
a  function  of  the  distance  r,  the  total  atmospheric  absorption  for  a  path  segment  of  length  d,  is 
formally  given  by: 

=  jS[H(r)]dr  (2-1'*) 

The  numerical  integration  of  Equation  2-14  is  performed  in  subroutine  ABSORB.  For  LOS 
paths  and  for  troposcatter  paths,  the  integration  is  by  Simpson’s  Va  rule  with  H(r)  and  therefore  7(H) 
determined  in  2-km  intervals  along  the  slant  path  r  (see  Figure  2-3).  For  a  diffracted  path,  the 
integration  in  Equation  2-14  for  the  two  paths  rj  and  Tj  within  the  horizons  of  the  transmitter  and 
receiver  is  obtained  in  the  same  manner.  In  the  diffracted  path  region  rj,  the  two  horizon  elevations 
hi  and  hj  are  averaged  and  the  integration  consists  simply  of  the  multiplication  of  the  value  7  for 
this  averaged  elevation  by  the  distance  between  the  two  horizons.  The  absorption  for  diffracted 
paths  is  then  the  sum  of  the  three  absorption  parts. 
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Troposcatter 


Figure  2-3.  Propagation  Paths  for  Different  Modes 
(Showing  H(r)  Determined  From  Path  Geometry) 
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2.5  TROPOSCATTER  PROPAGATION 

The  transmission  of  electromagnetic  energy  beyond  the  radio  horizon  occurs  by  tropospheric 
forward  scatter  as  well  as  by  diffraction.  In  this  section,  the  modeling  of  the  median  basic 
propagation  loss  resulting  from  tropospheric  forward  scatter  Ljro  is  discussed.  L^rq  is  given  by  an 
empirical  formula 


=  max  {[30  log  f  -  20  log  D  +  20  log  (r^/D)  +  F(D^)],  L,,} 

+  H  -  F„ 

0  0 


(2-15) 


where,  with  all  losses  in  dB 


Lps 

F(D0) 

Fo 

Ho 

To 

f 

D 

e 


free-space  loss 

attenuation  function 

scattering  efficiency  correction 

frequency  gain  function 

max  (D,D)  (see  Figure  2-4),  in  km 

frequency,  in  MHz 

great-circle  distance  between  transmitter  and  receiver,  in  km 
great-circle  arc  distance  between  transmitter  and  receiver,  in  radians. 


F(D6),  Fp,  and  are  empirical  functions  that  depend  on  angles  and  distances  derived  from 
the  terrain  intervening  between  transmitter  and  receiver.  These  functions  are  discussed  in 
Appendix  B. 


The  tropospheric  forward  scatter  loss  is  implemented  in  subroutine  TROPSC. 


^■^International  Radio  Consultative  Committee  (CCIR),  Tropospheric  Wave  Transmission  Loss 
Prediction.  CCIR,  Doc.  V/23-E,  March  1962. 
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In  Equations  2-17  and  2-18,  6^  ^  are  the  elevation  angles  (above  the  horizontal)  of  the  horizon 
obstacle  at  the  transmitter  (i  =  1)  and  at  the  receiver  (i  =  2).  These  are: 


tan  6, 


Vi 


K 


e,i 


dL.i 


L,i 


i  =  1,2 


(2-19) 


where 


hL,i  =  height  of  the  obstacle  that  produced  the  radio  horizon  above  MSL 

hs  i  =  structural  antenna  height,  plus  the  site  elevation  above  MSL 

dL,i  =  distance  from  the  antenna  to  the  radio  horizon 

a,  =  effective  earth  radius. 


Angular  elevation  of  horizon  rays  at  horizon  location  is: 


0,1 


(2-20) 


Horizon  obstacle  to  crossover  distances,  dcj,  are: 

dc,i  =  -  d^  j 

dc,2  =  ^OCi^^oA.2 


(2-21) 

(2-22) 


2.6  SPHERICAL-EARTH  LOSS 

Another  common  modeling  technique  used  in  both  TIREM  and  SEM  is  the  method  for 
calculating  and  combining  the  space- wave,  surface- wave,  and  diffracted- wave  losses  over  a  smooth 
spherical  earth.  This  process  is  programmed  in  subroutine  SEARTH  (and  its  ancillary  subroutines). 
Since  SEARTH  is  the  fundamental  basis  of  the  SEM  program,  its  theory  and  methodology  are 
described  in  the  next  section. 
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CHAPTER  3 

SPHERICAL-EARTH  MODEL  (SEM) 


3.1  GENERAL 

In  1941,  Norton^-'  published  a  systematic  method  for  calculating  the  ground- wave  field  for  both 
vertically  and  horizontally  polarized  waves  emitted  by  a  Hertzian  dipole  and  propagating  over  a 
spherical  earth.  Norton  summarized  the  results  of  a  number  of  classical  papers  dealing  with  the 
calculation  of  ground-wave  field  strength.  He  then  presented  a  method,  involving  a  combination  of 
simple  equations  and  graphs,  for  practical  calculation  of  ground- wave  field  strength.  The  method  is 
automated  for  use  in  TIREM  and  SEM  by  substituting  simple  equations  for  the  graphical  techniques 
of  Norton.  Essentially,  the  method  consists  of  calculating  the  space  wave  and  the  surface  wave, 
each  to  a  certain  distance,  and  the  diffracted  wave  beyond  a  critical  distance  by  the  first  term  of  the 
residue  series. By  making  a  plot  of  the  field  versus  the  distance,  the  intermediate  region  between 
the  two  distance  limits  can  then  be  traced  in  quite  easily.  The  formulation  given  here  and 
implemented  in  the  subroutine  SEARTH  (and  those  called  by  SEARTH)  is  virtually  the  same  as  that 
in  Norton’s  paper  with  the  following  exceptions: 

•  the  Hertzian  dipole  is  replaced  by  a  uniformly  radiating  source 

•  the  tracing  is  replaced  by  a  splice  function 

•  the  fields  are  referenced  to  a  free-space  field. 

This  chapter  provides  the  propagation  formulas  of  subroutine  SEARTH  in  four  sections  for  (1) 
the  space  wave,  (2)  the  surface  wave,  (3)  the  diffracted  wave,  and  (4)  the  splice  functions.  All  field 
formulas  are  normalized  to  the  free-space  field  Eo,  that  is,  the  field  identified  with  the  direct  wave. 


^■•Norton,  K.  A.,  "The  Calculation  of  Ground- Wave  Field  Intensity  Over  A  Finitely  Conducting 
Spherical  Earth,”  Proc.  I.R.E..  Vol.  29,  No.  12,  December  1941,  pp.  623-639. 

^■^Bremmer,  H.,  Terrestrial  Radio  Waves.  Elsevier  Publishing  Co.,  Amsterdam,  1949. 
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Before  presenting  the  formulas,  it  may  be  well  to  add  a  short  historical  note.  The  problem  of 
dipole  radiation  in  the  presence  of  a  dissipative  half-space  dates  back  to  Sommerfeld’s  classical 
paper^’^  and,  for  practical  purposes,  efforts  to  define  a  solution  continued  into  the  1960s.  There  are 
numerous  references  covering  the  subject.  As  a  starting  point,  the  interested  reader  should  consult 
the  books  authored  by  King,^"^  Bremmer  (Reference  3-2),  Banos, and  Wait.^'^  Bremmer  also 
presented  an  interesting  historical  survey  paper, and  in  the  same  publication  King  and  Wait^  * 
provided  a  consolidated  treatment  of  ground-wave  propagation. 

It  is  noted  here  that  the  propagation  formulas  in  ECAC’s  EFFSECC  model  (Reference  1-3)  are 
essentially  the  formulas  of  Reference  3-2.  In  the  diffracting  region,  a  great  number  of  the  terms  in 
the  residue  series  are  retained  and  consequently  the  propagation  subroutines  of  EFFSECC  require 
considerable  computer  time.  Because  of  a  requirement  for  a  fast  and  efficient  spherical-earth  model, 
the  SEM  model  was  implemented  based  on  Norton’s  paper  (Reference  3-1).  A  statistical 
comparison  between  SEM  and  EFFSECC  is  given  in  Chapter  5. 


^■^Sommerfeld,  A.,  "tfber  die  Ausbreitung  der  Wellen  in  der  Drahtlosen  Telegraphie,"  Ann.  Physics, 
28,  1909,  pp.  667-737. 

^■^King,  R.  W.  P.  and  Prasad,  S.,  Fundamental  Electromagnetic  Theory  and  Applications, 
Prentice-Hall,  Inc.,  Englewood  Cliffs,  NJ,  1986. 

^■^Banos,  A.,  Pinole  Radiation  in  the  Presence  of  a  Conducting  Half  Space.  Pergamon  Press,  Oxford, 
1966. 

^■^Wait,  J.  R.,  Electromagnetic  Waves  in  Stratified  Media.  Pergamon  Press,  Oxford,  1962. 

^■’Bremmer,  H.,  "A  Historical  Survey  of  the  Mathematical  Theories  of  Radio  Wave  Propagation," 
Symposia  Matematica  (Instituto  Nazionale  Di  Alta  Matematica,  Bologna,  Italy)  Vol.  18,  1976, 
pp.  57-83. 

^■®King,  R.  J.  and  Wait,  J.  R.,  "Electromagnetic  Ground  Wave  Propagation,  Theory  and  Experiment," 
Symposia  Matematica  (Instituto  Nazionale  Di  Alta  Matematica,  Bologna,  Italy)  Vol.  18,  1976, 
pp.  107-208. 
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3.2  SPACE-WAVE  PROPAGATION 

The  geometry  for  spherical-earth  calculations  at  points  within  LOS  is  shown  in  Figure  3-1. 
Constructive  or  destructive  interference  will  occur  between  the  direct  and  reflected  waves  at  the 
receiving  antenna  because  the  reflected  ray  path  is  electrically  longer  than  that  of  the  direct  ray. 

The  electrical  length  of  the  reflected  ray  path  is  greater  because  the  physical  path  length  is  longer 
and  because  of  the  phase  delay  that  occurs  at  the  point  of  reflection.  Also,  the  magnitude  of  the 
reflected  wave  will  be  smaller  than  that  of  the  direct  wave  because  it  is  proportional  to  the 
magnitude  of  the  reflection  coefficient  and  the  divergence  factor.  The  space  wave  is  formed  by  the 
vector  addition  of  the  direct  and  reflected  waves.  For  vertical  and  horizontal  polarization,  the 
space-wave  electric  field  Egpw,  normalized  to  the  direct  (or  free  space)  field  E„  at  the  field  point,  is 
expressed  by: 

lEspw^Ej  -  11  +  (r//r/)  D^R  eJ^l  D  <  d^p^  (3-1) 


where 


r/  =  +  (h/  -  h^O' 


(3-2) 


r/  =  +  (h/  +  h^O' 


(3-3) 


are  the  path  distances  traversed  by  the  direct  and  reflected  waves,  respectively. 


The  phase  factor  A  and  the  divergence  coefficient  D,,  are: 

A  -  (r2  rj ) 

A 


(3-4) 


2d,d2 

D,  =  [1  +  - ^ 


a  (hj  +  h2 ) 


(3-5) 


and  are  independent  of  polarization. 
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TRANSMlTTlHC  ANTENNA 


(FROM  REFERENCE  3-1) 

Figure  3-1.  Geometry  for  Spherical-Earth  Calculations  at  Points  Within  LOS 
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The  divergence  factor  is  the  ratio  of  field  strength  after  reflection  from  a  perfectly  reflecting 
spherical  surface  to  the  incident  field.  It  is  necessary  because  the  beam  reflected  from  a  spherical 
surface  diverges  at  a  greater  rate  than  it  would  if  reflected  from  a  plane  surface. 

The  heights  h/  and  hj'  are  the  antenna  heights  measured  from  the  tangent  plane; 

h,'  =  h,  -  (3-6) 

^  2  a. 


h^  - 


(D  - 

2  ae 


(3-7) 


where  hi  and  hj  are  the  structural  antenna  heights.  The  distance  di  to  the  reflection  point  is: 


di 


D 

2 


((l)/3  +  n/3) 


+  h^)  +  (D/2)2 


(3-8) 


which  is  the  trigonometric  solution  of  a  cubic  equation  with: 


cos  (J) 


3 

32a,  (hi  -  h^)D 


3 

4  [a,  (hi  +  h^)  +  (D/2)2]2 


(3-9) 
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The  reflection  coefficient  R  is  dependent  on  the  polarization  as  follows: 

For  Vertical  Polarization 


N^suh|;2^  -  -  cos^ijTj^ 

N^sunlfj^  +  -  cos^ijrj^ 


For  Horizontal  Polarization 


sini|r2'^  - 

sin\|r2^  +  yN^-co?ili^ 


with 


cos  t|;2^  =  D/r2'^ 

The  square  of  the  complex  index  of  refraction  is: 

+  j  60  oX 


where 


Cf  =  relative  dielectric  constant  of  earth 
a  =  earth  conductivity,  in  S/m. 


(3-10) 


(3-11) 


(3-12) 


(3-13) 
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The  normalized  space-wave  (geometric-optical)  field,  given  in  Equation  3-1,  oscillates  with 
distance  due  to  the  interference  of  the  direct  and  the  reflected  waves.  The  phenomenon  is  depicted 
in  Figure  3-2.  The  envelope  of  the  maximum  (C,  in  Figure  3-2)  is: 


and  the  envelope  of  the  minimum  (Cj  in  Figure  3-2)  is: 


Figure  3-2.  Numerical  Example  Showing  the  Interference  of  the  Direct  and  the  Reflected  Ray 
(f  =  100  MHz,  hi  =  hj  =  100  m,  e  =  10,  a  =  0.01  S/m,  Vertical  Polarization) 
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The  distance  to  the  last  maximum  of  the  oscillations  dwAx  the  distance  to  the  last  minimum  dMiN 
are  found  by  the  methods  in  Appendix  C.  The  distance  djp^  limiting  the  range  of  applicability  of 
Equations  3-1,  3-14,  and  3-15  lies  between  the  distance  of  the  last  maximum  d^Ax*  the  radio 
horizon  distance  dnoR,  where  the  geometric-optical  equations  fail,  and  is  given  by: 


^SPW  ^MAX  (^OR  *^MAx) 


(3-17) 


The  factor  Wj  =  0.6  is  empirically  determined  from  comparison  with  results  derived  with  the 
residue  series.  The  radio  horizon  distance  dnoR  is: 

^HOR  “ 

In  many  applications  it  is  preferable  to  replace  the  normalized  field  in  the  oscillating  region 
(Equation  3-1)  by  the  free-space  value  (Cq  in  Figure  3-2),  so  that  for  D  <  dgpw: 


l^sw/^ol  I^MAx/^oU 


(3-19) 


where  the  subscript  SPW  in  Egpw  has  been  replaced  by  FSW  to  distinguish  Equation  3-19  from 
Equation  3-1. 

When  the  free-space  wave  option  is  selected,  the  predicted  field  relative  to  free  space  is  given  by  the 
horizontal  line  Co  (see  Figure  3-2)  for  LOS  distances  in  the  oscillating  region.  At  distances  beyond 
the  intersection  of  Cq  and  C,,  the  field  is  given  by  the  curve  C^.  One  application  where  the  free- 
space  replacement  of  the  space  wave  is  required  is  in  TIREM.  The  use  of  a  geometric-optical  field 
based  on  smooth-earth  geometry  for  a  distance  less  than  d^Ax  would  be  unsuitable  in  a  rough-earth 
model  and  would  lead  to  large  prediction  errors. 
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3.3  SURFACE-WAVE  PROPAGATION 

The  ground  wave  consists  of  three  components:  the  direct  wave,  the  reflected  wave,  and  the 
surface  wave.  The  direct  and  reflected  waves  constitute  the  space  wave.  For  high  antennas  at  LOS 
distances,  the  direct  and  reflected  waves  may  combine  alternately  in-phase  or  out-of-phase  as 
distance  increases  causing  signal  peaks  and  nulls.  If  the  distance  continues  to  increase  until  the 
receiving  antenna  approaches  the  horizon,  the  lengths  of  the  direct  and  reflected  paths  become 
nearly  the  same,  and  because  of  a  phase  reversal  of  the  reflected  wave,  the  direct  and  reflected 
waves  nearly  cancel  leaving  only  the  surface  wave.  The  surface  wave  is  of  primary  importance 
when  both  antennas  are  within  a  few  wavelengths  of  the  ground.  For  low  antennas,  less  than  about 
a  wavelength,  the  surface  wave  can  be  significant  compared  to  the  space  wave  at  short  LOS 
distances.  The  surface-wave  field  intensity  normalized  to  the  free-space  intensity  is: 


E 


SUW' 


/Eo 


'  |(l-R)l(r//r20  f(p)  P<0.729 

|(1-R)|  min  f(p)5  pkO.729 


where  f(p)  is  the  Sommerfeld  reduction  factor  and  Fj  ^  are  functions  described  below.  The  ratio 
f(p)/p  is  plotted  in  Figure  2  of  Reference  3-1.  The  Sommerfeld  reduction  factor  f(p)  is 
approximated  here  by  the  empirical  relation: 


f(p) 


1  +  3p 

1  +  3p  +  6p^ 


(3-21) 


The  parameter  p  is  the  Sommerfeld  numerical  distance  and  is: 


TZ  +  (h^  +  h^)" 


(3-22) 


where,  y  =  Vi  for  vertical  polarization  and  y  =  -Vi  for  horizontal  polarization. 
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The  function 


■SUW,i 


1  + 

2nh.  'f 

N 

[aIN^k] 

(3-23) 


is  a  height-gain  function  for  a  surface  wave  area  over  a  plane  earth. 

The  distance  djuw,  limiting  the  range  of  applicability  of  the  surface  wave  formulas,  is  generally 
taken  to  be: 


dsuw  =  12000  A’ 

3.4  DIFFRACTED-WAVE  PROPAGATION 

For  low  antennas  at  a  short  distance,  the  direct  and  ground  reflected  waves  cancel  leaving  only 
the  surface  wave.  The  surface  wave  at  large  distances  is  a  diffracted  wave.  The  curvature  of  the 
earth  prevents  the  waves  from  passing  directly  into  a  shadow  region.  They  must  reach  this  shadow 
region  by  a  process  of  bending  around  the  curved  surface  of  the  earth.  This  process  is  called 
diffraction. 

For  distances  D  >  doiw,  the  diffracted  field  intensity  Fdiw  normalized  to  the  free-space  field  is 
given  as  the  first  term  of  the  residue  series: 

lEon^J  =  113  Fj  F,  I  VV  exp[-(2jt)'«il'l  (3-25) 


where 
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r\'  =  PtiD 
ri  = 


(3-26) 

(3-27) 


Y  =  min  (0.03,  0.08  K^} 


7  is  an  approximation  to  the  plot  of  Figure  6  of  Reference  3-1. 


The  height-gain  functions  Fj  in  Equation  3-25  (for  i  —  1,2)  are: 


F.  = 


max  (Fgy^  j,  Fjjj^^}  ,  hj^l 

^SUW,i  ’ 


The  height-gain  function  for  elevated  antermas  BLOS  is: 


where  fj  (hj)  is  a  polynomial  approximation  to  Figure  11  in  Reference  3-1. 


The  parameters  hj  and  8  are  defined  as: 


h:  = 


P" 


(a.AO 


2n1/3  * 


6  =  max  (0.328776/K,  0.357945) 


(3-28) 


(3-29) 


(3-30) 


(3-31) 


(3-32) 
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The  parameter  /3  in  Equations  3-25  and  3-26  is: 


p  =  0.9[e-^26K2  +  0.78] 


(3-33) 


where 


K  =  (27iayX)-^/3|N2|y 


(3-34) 


Foiw.i  is  Equation  20  of  Reference  3-1  and  8  and  /3  are  approximations  to  the  plots  of  Figures  10 
and  5,  respectively  of  Reference  3-1.  The  distance  dow  is: 


^DIW  ^HOR 


1.5 

Pq 


(3-35) 


3.5  SPLICE  FUNCTIONS 

For  situations  where  the  minimum  distance  doj^v,  beyond  which  the  diffraction  calculation 
applies,  is  greater  than  the  maximum  distances  djpw  or  dsu^,  where  either  the  space  wave  or  surface 
wave  calculations  apply,  it  is  necessary  to  estimate  the  normalized  field  by  using  a  splice  function. 
This  function  is  equivalent  to  Norton’s  graphical  method  (Reference  3-1,  p.  632)  in  which  he 
suggests  that  "a  smooth  transition  curve  may  be  drawn. "  The  splice  function  |  Espuce/Eo  |  developed 
for  this  purpose  simulates  the  normalized  field  between  these  limits. 

The  technique  is  illustrated  in  Figure  3-3.  The  figure  is  divided  into  three  regions.  In  the 
leftmost  region,  the  normalized  field  is  represented  by  |  \  (D)  and  can  be  either  |  Egpw/Eo  |  or 

I  Esuw/Eo  I  depending  on  which  is  dominant.  This  region  ends  at  a  distance  of  either  dspw  or  djuw. 
given  generically  as  Dj.  The  rightmost  region  is  the  diffraction  region  where  the  field  is  represented 
by  IEdpw/EoKD).  This  region  begins  at  Do  =  dew  The  splice  function  in  the  middle  region  serves 
to  connect  |Es/Eo|(D)  to  |Ediw/Eo|(D). 
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Figure  3-3.  Example  of  the  Splice  Function  Joining  Either  the  Space-Wave  or  Surface-Wave  Fields 

with  the  Diffracted- Wave  Field 


The  splice  function  is  of  the  form 


IEspucb'E.KD)  =  |E^J(D,)  e'-®’’-®’' 

where  lEs/EjCDs),  D^,  and  are  determined  such  that  the  normalized  field  is  continuous  at  the  pair 
of  limiting  distances  Dj  =  dgpw  or  dju^,  as  appropriate,  and  Du  =  doiw 

The  parameter  jSs  is; 


Ps  = 


1 


■*DIW 


|Es/EJ(D,) 

|Ed™/E.I(Dd) 


(3-37) 
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Equation  3-36  is  the  general  form  for  the  splice  function.  The  coefficient  |  Ej/Eq  |  (D^)  is  either 
the  normalized  space-wave  field  at  =  dsp^  or  the  normalized  surface  wave  field  at  =  dju^,  as 
appropriate.  In  Equation  3-37,  |  Ed,w/Eo  |  (Do)  is  the  normalized  diffracted  wave  field  evaluated  at 
the  limiting  distance  =  dow  The  splice  function  makes  the  field  strength  continuous  at  the 
transition  points.  However,  the  first  derivative  of  the  field  strength  is  discontinuous  at  the  transition 
points. 

3.6  SUMMARY  OF  SPHERICAL-EARTH  PROPAGATION  FORMULAS 


The  excess  propagation  loss  is  calculated  by  multiplying  the  logarithm  of  the  normalized  field  by 
-20.  The  excess  propagation  loss  for  the  space  wave,  free  space,  surface  wave,  splice  function,  and 
diffracted  waves  are  designated  as  Aspw,  Aps^,  Agu^^  Aspuce>  respectively.  Table  3-1 

contains  a  summary  of  the  excess  loss  calculations,  the  range  of  applicability,  and  a  cross  reference 
to  the  equation  in  the  text.  The  SEM  excess  loss  Asem  is  the  excess  loss  for  the  appropriate  wave 
type  depending  on  distance.  At  distances  less  than  dow,  it  is  possible  that  the  formulations  for  both 
the  space  wave  and  the  surface  wave  are  applicable.  In  such  a  case,  the  appropriate  loss  is  the 
smaller  of  Agpw  or  Ajuw  The  total  spherical-earth  propagation  loss  calculated  by  SEM  is: 


I  Lps  ^SEM  Aabsorb  paths 

I  min  [Lpg  +  Ag^j^  +  Ay^gQj^g,  L.j.j^q]  Beyond  LOS  paths 


TIREM  also  requires  the  calculation  of  a  spherical-earth  propagation  loss  relative  to  free  space 
Aseloss'  Aseloss  is  the  appropriate  loss  of  the  following  subset  from  Table  3-1:  Apsw.  Ajuw, 

Asplice’  ^nd  Adiw  Use  of  the  theoretical  space  wave  Agpw.  calculated  for  a  smooth  earth,  or  its 
upper  or  lower  envelopes,  Amax  or  Amin,  would  be  inappropriate  for  a  rough-earth  model  and  would 
contribute  to  excessive  prediction  errors  for  LOS  paths.  These  optional  modes  are  therefore  not 
used  in  TIREM. 
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Table  3-1.  Summary  of  Spherical-Earth  Propagation  Formulas 


Wave  Type 

Excess  Propagation  Loss 

Region  of 
Applicability 

Equation  in 
Text 

Space  Wave 

^SPW  ~  1  ^spy/ZEj,  1 

D  <  dspw 

3-1 

Envelope  of  Space  Wave 

^MAX  "  1  EmaX^Eo  I 

^MIN  ~  1  I 

D  ^  d^iAx 
^  dMiN 

3-14 

3-15 

Space  Wave 

Approximation 

D  <  dgpyr 

3-19 

Surface  Wave 

D  ^  dsuw 

3-20 

Space/Diffracted 

Wave  Splice 

Asplice  -  -20  log  1  Espuce/Eo  1 

^SPW  ^  D  ^  ^DIW 

3-36 

Surface/Diffracted 

Wave  Splice 

^suw  D  <  df)!^ 

Diffracted  Wave 

AdIW  “  i  Ediw/Eo  I 

D  >  dppw 

3-25 

The  equations  defining  the  distances  are  as  follows:  djpw.  (Equation  3-17);  djuw.  (Equation  3-24);  dow, 

(Equation  3-35). 
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CHAPTER  4 

TERRAIN-INTEGRATED  ROUGH-EARTH 
MODEL  (TIREM) 


4.1  GENERAL 

The  Terrain-Integrated  Rough-Earth  Model  (TIREM)  is  designed  for  calculation  of  the  reference 
basic  median  propagation  loss  (path  loss)  using  digital  terrain-elevation  data  (DTED)  to  construct  the 
terrain  profile  along  the  great-circle  path  between  two  antennas.  A  combination  of  theoretical 
techniques  and  empirically-derived  equations  is  used  in  TIREM  to  produce  acceptable  estimates  of 
path  loss  over  irregular  terrain  for  most  purposes.  These  methods  are  integrated  formally  into  a 
TIREM  module,  which  has  been  validated  by  comparison  with  measured  propagation  data. 

Statistical  comparisons  of  predicted  propagation  losses  with  corresponding  measured  losses  are 
provided  in  Chapter  5. 

4.2  DIGITAL  TERRAIN  DATA  IN  PROPAGATION  MODELING 

The  modeling  of  propagation  in  TIREM  depends  on  terrain-elevation  information.  TIREM 
requires  a  terrain  profile  described  by  a  set  of  discrete  points,  the  position  of  which  is  specified  by  a 
distance  from  the  transmitter  and  an  elevation  above  mean  sea  level.  TIREM  evaluates  each  such 
profile  to  determine  whether  propagation  takes  place  over  land  or  sea  or  a  combination  of  the  two. 
TIREM  differentiates  land  from  sea  by  noting  that  profile  elevations  for  land  are  not  equal  to  0 
while  elevations  for  sea  are  equal  to  0.  The  user-specified  earth  conductivity  and  relative 
permittivity  are  used  if  the  path  is  over  land  and  the  corresponding  characteristics  for  sea  water  are 
used  if  the  path  is  over  sea.  These  electrical  characteristics  for  land  and  sea  differ  significantly  and 
greatly  influence  the  predicted  propagation  loss  for  smooth  or  nearly  smooth  surfaces. 
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Using  effective  earth  radius  geometry  and  the  terrain  profile,  TIREM  then  determines  whether 
the  antenna  terminals  are  within  LOS  or  BLOS.  This  is  done  by  calculating  the  slope  (tangent  of 
the  elevation  angle)  to  each  profile  point  as  measured  at  the  transmitting  antenna.  The  calculation  is 
done  in  subroutine  SLTANG,  using  Equation  2-19. 

The  process  begins  by  assuming  that  the  maximum  slope  is  that  of  the  receiving  antenna  (i.e., 
that  the  path  is  LOS).  Then,  the  slope  to  each  profile  point  between  the  antennas  is  calculated 
beginning  at  the  receiver  and  working  back  to  the  transmitter.  If  the  slope  to  the  receiving  antenna 
is  greater  than  that  to  all  other  profile  points,  the  path  is  LOS,  otherwise  it  is  BLOS. 

If  the  path  is  BLOS,  the  profile  point  forming  the  steepest  slope  is  recognized  as  the  transmitter 
horizon.  That  point  then  becomes  a  beginning  point  of  observation  for  the  measurement  of  slopes  to 
subsequent  profile  points.  If  that  observation  point  is  within  LOS  of  the  receiving  antenna,  the 
transmitter  and  receiver  have  a  common  horizon.  If  not,  subsequent  points  with  the  highest 
elevation  angles  become  successive  points  of  observation  (and  are  also  identified  as  diffraction  kmfe 
edges)  until  an  observation  point  is  found  that  is  within  LOS  of  the  receiving  antenna.  That  point  is 
then  the  receiver  horizon. 

In  the  majority  of  cases,  the  required  terrain  information  is  taken  from  extensive  coverage  of 
level  1  digital  terrain-elevation  data  generated  by  the  Defense  Mapping  Agency  (DMA).'*  *  DMA 
provides  Digital  Terrain-Elevation  Data  (DTED).  The  DTED  level  1  format  is  a  uniform  matrix  of 
terrain-elevation  values  spaced  at  3  arc-second  (approximately  100  m)  intervals  to  approximate  the 
vertical  resolution  of  1:250,000  scale  elevation  contour  maps.  The  terrain-elevation  values  are  not 
exact.  The  linear  error  in  the  elevation  of  a  point  with  respect  to  mean  sea  level  is  ±30  m;  the 
circular  error  of  a  point  in  the  horizontal  position  with  respect  to  the  World  Geodetic  System  is  50 
m  (both  errors  are  at  90-percent  confidence  level). 


‘*'^Militarv  Specification  For  Digital  Terrain  Elevation  Data  (DTED)  Level  1  (ICD),  MIL-D-89000, 
Defense  Mapping  Agency,  Fairfax,  VA,  26  February  1990. 
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Studies  conducted  at  EC  AC  established  that  for  the  purpose  of  propagation  modeling,  thinning 
the  data  by  a  factor  of  5  to  a  resolution  of  15  arc  seconds  did  not  degrade  any  of  the  performance 
measurements  considered/'^  With  this  15  arc-second  grid  spacing,  the  storage  requirement  for  the 
Level  1  DIED  is  reduced  by  a  factor  of  25  over  that  required  for  archiving  the  3  arc-second 
resolution  data  provided  by  DMA.  A  further  reduction  of  storage  space  by  a  factor  of  about  3:1  is 
effected  by  quantization  of  the  elevation  into  3-meter  (10-foot)  increments  rather  than  the  1 -meter 
increments  of  the  raw  Level  1  DIED.  ECAC  studies  have  shown  that  using  the  3-meter  elevation 
quantum  does  not  degrade  TIREM  error  statistics  when  comparing  model  predictions  with 
measurements.  The  ECAC  operational  topographic  data  system  therefore  uses  a  3-meter  elevation 
quantum  and  a  15  arc-second  spatial  resolution  for  archiving  the  Level  1  DIED. 

4.3  PROPAGATION  OVER  LAND 

4.3.1  I.ine-of-Sight  Region 


For  propagation  over  irregular  earth  in  the  LOS  region,  there  is  coupling  between  the  antennas 
because  of  the  direct  free-space  wave.  Also,  there  is  a  coupling  path  via  an  earth-reflected  wave; 
however,  this  path  may  or  may  not  be  important.  Unlike  the  LOS  case  discussed  in  the  previous 
section  for  propagation  over  a  smooth  spherical  earth,  there  generally  is  no  broad  expanse  of  level 
terrain  to  provide  a  single  significant  reflected  wave  at  the  receiver.  Instead,  the  earth  resembles  a 
multi-faceted  solid  composed  of  reflecting  planes  of  many  sizes,  shapes,  and  orientations.  It  is,  of 
course,  impractical  to  represent  such  an  earth  by  a  terrain  profile  confined  to  the  great-circle  plane, 
as  TIREM  attempts  to  do.  For  this  reason,  at  frequencies  above  20  MHz,  TIREM  does  not  attempt 
to  predict  ground- wave  propagation  in  the  manner  of  SEM.  Instead,  TIREM  predicts  a  free-space 
loss  unless  it  determines  that  terrain  protrudes  into  the  first  Fresnel  zone. 


'‘-^Cameron,  S.  and  Kuebler,  W.,  "Digital  Terrain  Data  in  Radio  Wave  Propagation  Predictions: 

An  Assessment,"  Proc.  DMA  Symposium  89,  DMA  Systems  Center,  Reston,  VA,  15-17  May  1989 
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The  first  Fresnel  zone  is  defined  as  the  spheroid  of  revolution  about  the  line  joining  the  antennas 
(the  direct  path)  such  that  the  sum  of  the  distances  from  any  point  on  its  surface  to  each  of  the 
antennas  exceeds  the  direct  path  by  exactly  one-half  wavelength.  If  sufficient  penetration  into  the 
Fresnel  zone  region  occurs,  an  additional  space-wave  reflection  loss  is  added  to  the  free-space  loss. 

Because  of  the  requirement  that  the  reflection  point  be  within  the  first  Fresnel  zone,  it  is  possible 
to  make  small  grazing  angle  approximations  without  loss  of  validity.  It  is  also  possible  to  assume 
that  the  antenna  gains  associated  with  the  direct  and  reflected  field  contributions  within  the  first 
Fresnel  zone  are  virtually  the  same.  It  is  this  concept  that  makes  possible  the  calculation  of  the  so- 
called  basic  propagation  loss  between  isotropic  antennas. 

In  the  LOS  region,  the  terrain  intercepts  a  portion  of  the  radiation  emitted  by  a  source,  and  part 
of  the  intercepted  energy  is  reflected  in  a  direction  leading  to  a  receiving  point  as  shown  in  Figure 
4-1.  Interference  then  takes  place  between  the  radiation  arriving  via  the  direct  ray  and  that  arriving 
via  the  reflected  ray.  The  space  wave  Fgpw  normalized  to  the  direct  wave  F^  for  an  irregular  surface 
is  given  by: 

lEs^/EJ  = 


This  equation  is  the  same  as  Equation  3-1  for  the  special  case  of  rj’  «  rj’  (when  both  antennas 
are  near  the  ground)  except  that  the  divergence  coefficient  is  replaced  by  the  scattering  coefficient 
p„.  This  equation  applies  for  reflection  from  a  rough  planar  surface  that  approximates  the  dominant 
reflecting  surface  in  irregular  terrain.'^'^ 


'‘■^Schelleng,  J.,  Burrows,  C.,  and  Ferrell,  E.,  "Ultra-Short- Wave  Propagation,"  Proc.  I.R.E..  Vol.  21, 
No.  3.,  March  1933,  pp  427-463. 
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d 


(EXAMPLE  A)  INTERCEPTED  BY  A  FLAT  TERRAIN  PLANE 


d 


(EXAMPLE  B)  INTERCEPTED  BY  UNDULATING  TERRAIN 

Figure  4-1.  Geometry  of  the  First  Fresnel  Zone  Spheroid 
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In  the  modeling,  Equation  4-1  is  approximated  by  a  first-order  expansion  of  cos  A,  so  that: 


iEspw/Ej  =  /l  +  -  2  cos  A 


(4-2a) 


^(1  -  A^ 


(4-2b) 


These  formulas  are  based  on  the  ray  representation  implying  that  the  principal  part  of  the 
specularly  reflected  energy  arises  from  an  area  near  the  point  of  reflection  of  the  maximum  terrain 
intrusion  into  the  first  Fresnel  zone.  The  total  reflection  area,  in  fact,  extends  to  an  area  of 
approximately  half  the  area  formed  by  the  intersection  of  the  first  Fresnel  zone  spheroid  with  the 
terrain.  Equations  4-2  are  correct  only  for  flat  terrain.  For  undulating  terrain,  they  must  be 
considered  as  approximations.  Note,  however,  that  the  space-wave  formulation  does  not  depend  on 
the  polarization  of  the  source  nor  on  the  dielectric  constant  or  conductivity  of  the  terrain.  The 
approximation  given  in  Equation  4-2b  is  actually  used  in  TIREM. 

The  phase  factor  A  for  terrain  points  i  =  1,  2,  ...  is,  as  shown  in  Appendix  D,  expressed  in 
terms  of  the  clearance  C;  and  the  width  of  the  first  Fresnel  zone  f;  as: 


A  =  IT  min  {(c/fi)^  1} 


(4-3a) 


*  X  min  {(c/fi)^  -} 

TT 


(4-3b) 
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The  approximation  given  in  Equation  4-3b  corresponds  to  the  approximation  of  Equation  4-2b 
and  is  used  in  TIREM.  The  clearance  Cj  is  the  distance  of  the  direct  ray  above  the  i-th 
terrain-elevation  point  (see  Figure  4-1)  and  fj  is  given  by: 


fi  = 


x/(^ " 


(4-4) 


where  dj  i  and  ^2  i  are  the  direct  ray  distances  from  the  i-th  terrain  point  to  the  transmitter  and 
receiver,  respectively.  Use  of  the  factor  I/tt  in  Equation  4-3b  restricts  the  phase  factor  to  A  <  1  so 
that  Espw/Eo  <  1  and  Equation  4-2b  applies  only  if  (cM  ^  For  larger  clearance  ratios,  the 
substitution  A  =  1  so  that  Ejpw/E,,  =  1  is  made,  thereby  replacing  the  oscillating  lobe  structure  of 
the  normalized  field.  Equation  4-1,  by  the  free-space  value.  Refer  to  Appendix  D  and  Figure  D-2 
for  a  complete  explanation  and  an  illustration  of  this  substitution.  Although  the  oscillating  structure 
of  the  field  is  observed  in  practice  when  the  terrain  is  relatively  flat,  modeling  of  the  field 
oscillations  is  not  advisable.  The  predicted  oscillations  are  nearly  always  shifted  with  respect  to  the 
observed  field  and  thus  result  in  large  modeling  errors  despite  the  fact  that  the  oscillating 
characteristic  of  the  field  is  preserved.  Consequently,  it  proved  both  practical  and  expedient  in 
modeling  to  eliminate  the  oscillation.  Equation  4-3b  also  implies  that  the  clearance 
C.  =  yfi^  ^  0.56  fi  (the  terrain  effects  on  the  propagating  wave  are  restricted  to  clearances 

well  within  the  first  Fresnel  zone). 


The  reflection  coefficient  R  (given  by  Equations  3-10  and  3-11)  is  very  nearly  -1  regardless  of 
polarization  when  the  angle  between  the  reflected  ray  and  the  ground  is  small,  such  as  is  the  case 
when  the  reflecting  terrain  points  are  well  within  the  first  Fresnel  zone.  The  scattering  coefficient 
Po  represents  the  effects  of  the  irregularities  of  the  ground,  and  a  mean  value  has  been  established 
empirically  from  measured  propagation  data  and  depends  upon  the  wavelength  X  as: 


Po 


=  e 


-xyx 


(4-5) 


with  \  =■  0.02  meters. 
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The  loss  relative  to  free  space  due  to  the  reflected  ray  Aref  is  evaluated  in  subroutine  RFLECT 
and  is; 


Aref  =  -20  log(  |  E^p^/Ej ) 

=  -10  log  [(1  -  +  p„  A^] 


(4-6) 


At  grazing  incidence,  where  the  phase  factor  A  is  zero,  Equation  4-6  is  expressed  as: 


A„  =  Apgp  (A  =  0)  =  -20  log  (1  -  p„) 


(4-7) 


The  above  formulation  accounts  for  the  space-wave  part  of  the  LOS  field  contribution.  For 
frequencies  below  approximately  20  MHz  where  the  wavelength  is  15  meters  or  more,  moderate 
terrain  irregularities  become  less  important  and  it  is  reasonable  to  approximate  LOS  propagation 
using  the  space-  and  surface-wave  formulations  described  in  the  previous  section  for  propagation 
over  a  spherical  earth.  This  is  done  in  TIREM  by  employing  the  spherical-earth  algorithms 
exclusively  for  frequencies  below  16  MHz  and  by  combining  the  spherical-earth  loss  and  reflection 
loss  by  linear  interpolation  for  frequencies  between  16  and  20  MHz. 

The  surface-wave  field  contribution  is  estimated  through  the  field  calculation  based  on  a 
spherical  earth  (no  simple  method  exists  to  derive  the  surface-wave  field  in  the  presence  of  terrain 
features)  and  is  given  by  Equation  3-20.  The  surface-wave  calculation  is  part  of  subroutine 
SEARTH  (refer  to  paragraph  3.5  and  Table  3-1),  resulting  in  a  loss  relative  to  free  space  Aseloss. 
and  depends  on  antenna  heights,  distance,  and  dielectric  constant  and  conductivity  of  the  terrain. 
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The  LOS  loss  relative  to  free  space  is  then  estimated  to  be: 


A 


LOS 


min  (Ag^Loss’  ^ref) 

f 

^MHz 

< 

16 

min  [a  A^g^oss  (1“^)  Aj^p,  Aj^p] 

16  ^  fMH. 

< 

20 

A 

^REF 

ImHz 

> 

20 

(4-8) 

where 


(4-9) 


4.3.2  Diffractinp  Region 

In  the  region  where  the  radiation  emitted  by  a  source  is  diffracted  by  prominent  terrain  features 
such  as  hills,  these  features  are  modeled  as  ideal  knife-edges.  Reflection  may  also  be  present. 
Figure  4-2  shows  the  situation  where  the  field  is  diffracted  by  k  ridges.  In  this  figure,  the 
terminating  points  (transmitting  and  receiving  antennas)  are  represented  by  darkened  circles  and  the 
diffracting  points  by  open  circles.  For  k  knife  edges  the  field  intensity  Edif  normalized  to  the  free- 
space  intensity  is: 

IWEol  =  «l.K(Pk)«2.k  (4-10) 


k=l 


Figure  4-2.  Diffraction  Over  Two  Knife  Edges,  Epstein-Peterson  Construction 
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The  knife-edge  diffraction  term  F,;  for  each  knife  edge  is  the  magnitude  of  the  Fresnel  integral: 


Aexp  (“V^u^  +  0.252Uk) 
1 

27m, 


0  <  u,  <  1.8 
u,  >  1.8 


(4-11) 


The  Fresnel  integral  in  Equation  4-11  is  approximated  in  subroutine  SKEDIF  by  a  convenient 
curve  fit  for  values  of  U;.  between  0  and  1.8  and  by  its  asymptotic  value  for  u^  >  1.8.  The 
dimensionless  parameter  u^  is  given  by: 


Uk  =  h,/f, 


(4-12) 


where  h^  is  the  "height"  of  the  k-th  knife  edge  as  shown  in  Figure  4-2  and: 


(4-13) 


is  the  radius  of  the  first  Fresnel  spheroid  at  the  k-th  knife  edge,  and  dA.k  and  dg.k  are  the  separation 
distances.  The  parameter  u^  is  related  to  the  frequently  used  parameter  v  by  v  =  Vl  %  (Reference 
1-4,  p.  7-1). 

For  a  single  ridge  where  k  =  1,  Equation  4-10  reduces  to  an  approximation  for  the  four-ray 
expression  given  first  by  Bullington.'*'^  The  diffracting  ridge  is  treated  as  a  kmfe  edge,  and 
reflection  can  occur  in  the  reflecting  regions  between  the  transmitting  point  and  the  ridge  and 


'’Radio  Wave  Propagation:  Consolidated  Summary  Technical  Report  of  the  Committee  on  Propagation 
of  the  National  Defense  Research  Committee.  Academic  Press,  New  York,  NY,  1949,  p.  69. 
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between  the  ridge  and  the  receiving  point/'^  For  more  than  one  ridge,  Equation  4-10  incorporates 
the  extended  Epstein-Peterson  method'*'^  wherein  the  first  ridge  is  treated  as  the  source  illuminating 
the  second  ridge,  etc.  Terrain  effects  between  the  ridges  are  not  taken  into  account;  hence,  only  the 
region  preceding  the  first  peak  and  the  region  following  the  last  peak  are  considered  reflecting 
regions.  The  contributions  ai  ^  and  a2,k  reflecting  regions  result  from  the  phase  shifts  in  the 

rays  due  to  specular  reflection  and  are  expressed  as  in  the  LOS  case.  Equation  4-2b,  by: 


«i.k  =  \/(l  -  P,f  +  PkA^ 


(4-14) 


«2,k  =  \/(l  -  Pk)"  +  pA 


(4-15) 


The  phase  factors  A)  and  Aj  for  the  two  reflecting  regions  are  evaluated  as  in  Equation  4-3b. 
The  mean  scattering  coefficient  is  empirically  determined  from  measured  propagation  data  and 
has  the  same  wavelength  dependence  as  in  the  LOS  case.  The  coefficient  p^  reduces  to  the  LOS 
case.  Equation  4-5,  for  k  =  0,  and  is  given  by; 


^k 


g-(ktl)\„/X 


(4-16) 


From  an  examination  of  the  terrain  profiles  where  diffraction  takes  place,  and  associated 
measured  propagation  data,  it  appears  that  the  reflecting  regions  are  not  on,  or  near,  the  slopes 
leading  to  the  point  of  diffraction.  Hence,  those  terrain  points  within  5  percent  of  the  distance  from 
the  ridge  to  the  terminal  point  are  excluded  from  the  calculation  of  Aj  and  A2.  It  is  also  evident  that 


'‘■^Crysdale,  J.  H.,  "Large  Reduction  of  VHF  Transmission  Loss  and  Fading  by  the  Presence  of  a 
Mountain  Obstacle  in  Beyond-Line-of-Sight  Paths,"  Proc.  IRE,  Vol.  43,  No.  5,  May  1955,  p.  627. 

'•■^Epstein,  J.  and  Peterson,  D.  W.,  "An  Experimental  Study  of  Wave  Propagation  at  850  MC," 

Proc.  IRE.  Vol.  41,  No.  5,  May  1953,  pp.  595-611. 
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such  a  restriction  is  required.  Otherwise,  an,  and  a2,i,  would  always  be  evaluated  by  their  limiting 
value  A,  =  A2  =  0  because  points  close  to  the  ridge  obstruct  the  clearance. 

Extensive  evaluation  of  predictions  using  Equation  4-10  versus  measured  propagation  data 
showed  that  the  best  results  were  obtained  by  limiting  the  specular  reflection  to  the  one  region  for 
which  the  effect  of  the  phase  relations  between  an,,  a2j„  and  is  to  reduce  the  received  field 
intensity  the  most.  In  the  modeling,  therefore.  Equation  4-10  is  replaced  by: 


(4-17) 


It  is  noted  that  the  above  knife-edge  diffraction  equation  is  independent  of  the  conductivity  and 
dielectric  constant  of  the  terrain.  The  equation  is  applicable  for  frequencies  greater  than  20  MHz 
and  is  evaluated  in  subroutine  DIF.  The  result  is  the  knife-edge  diffraction  loss  above  free  space 

-^KNIFE- 


•'^KNIFE  ^°§(  i  ^DIF'^^o  I  ) 

=  max(Aj  „A2,)  +  Ap, 


(4-18) 


where 


Ai,k  =  -20  log  a, , 


(4-19) 


A2,k  =  -20  log  a^ , 


(4-20) 


Ap ,  =  -20  log  F, 


(4-21) 
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At  grazing  incidence,  where  the  phase  factors  Aj  or  Aj  are  zero,  Equations  4-19  and  4-20  are 
identical,  and  using  Equations  4-14  and  4-15  may  be  expressed  as: 

A^  =  0!j  ^  (Aj  =  0)  =  (^2  ^  ^  log  (1  ■  Pk)  (4-22) 


As  in  the  LOS  case,  the  exeess  loss  defined  above  for  grazing  incidence  A^  is  required  in 
paragraph  4.3  for  land/sea  propagation.  Moreover  for  k  =  0,  this  definition  agrees  with  Equation 
4-7. 


The  diffraction  loss  is  as  given  by  Aknife  except  for  terrain  features  approaching  the  smoothness 
of  a  spherical  earth.  The  modeling  in  such  a  case  would  detect  a  knife  edge  at  every  terrain  point 
resulting  in  unreasonable  high-loss  values,  as  every  diffraction  point  would  result  in  a  loss  of  at  least 
6  dB.  Consequently,  if  there  are  more  than  two  knife  edges,  an  average  knife-edge  loss  value  of 
Ake  is  calculated  by: 


^KE  ^^F,k 

The  diffraction  loss  for  20  MHz  and  above  is  given  by: 


(4-23) 


A 


DIF.VHF 


min(A5gLoss’''^KNiFE)  k>3  and  A^g  <  7dB 
Aknife  k<3  or  \g  >  7dB 


(4-24) 


where  A^e  =  7  dB  is  an  empirically  established  transition  value.  Aseloss  Is  the  spherical-earth 
estimate  of  the  surface-wave  and  diffracted-wave  contribution  and  depends  on  the  relative 
permittivity  and  conductivity  of  the  terrain.  Aseloss  is  evaluated  in  subroutine  SEARTH  (refer  to 
paragraph  3.5  and  Table  3-1).  If  there  are  fewer  than  three  knife  edges,  the  diffraction  loss  is 

simply  Aknife- 
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For  frequencies  less  than  16  MHz,  the  diffraction  loss  can  be  estimated  by  the  spherical-earth 
loss,  as  shown  in  the  statistics  in  Chapter  5.  Hence,  for  frequencies  from  1  through  16  MHz: 


A 


DIF.HF 


‘SELOSS 


W  ^  16 


(4-25) 


The  loss  for  the  intermediate  frequency  region  16  <  f^Hz  ^  20  is  calculated  by  a  linear 
interpolation: 


A 


DIF, a 


^  ■'^DIF.HF  (1  ■'^DIF.VHF 


(4-26) 


where  a  is  given  by  Equation  4-9. 

4.3.3  Summary  of  Land  Propagation  Formulas 

The  path  loss  predicted  by  TIREM  consists  of  the  appropriate  loss  terms  defined  in  Chapter  2 
and  paragraph  4.3  of  this  handbook  and  depends  on  whether  the  path  is  LOS  or  BLOS.  For  a 
BLOS  path,  the  total  loss  for  both  diffraction  and  troposcatter  are  calculated  and  the  dominant 
coupling  (minimum  of  the  two  path  losses)  is  assumed  to  prevail.  The  total  terrain-dependent 
propagation  loss  for  LOS  and  BLOS  paths  is: 


4s  *  *  A„^„  LOS  (4-27) 

min  (Lp5  +  Apjp  +  A^g^Q^g,  Lj-rq  ■'^absorb)  BLOS 


The  excess  propagation  loss  as  obtained  from  the  normalized  fields  for  LOS  or  diffraction  is 
shown  in  Table  4-1.  Alqs  is  the  appropriate  excess  loss  for  line-of-sight  as  determined  by  the 
frequency.  A^jp  is  the  appropriate  excess  loss  for  the  diffracting  region.  Absorb  is  the  molecular 
absorption  loss  calculated  by  Equation  2-14  as  appropriate  for  each  case  shown  in  Figure  2-3.  The 
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LOS,  knife-edge-diffracted,  and  spherical-earth-diffracted  propagation  losses  are  calculated  in 
subroutines  TIRLOS,  DIF,  and  SEARTH,  respectively.  The  troposcatter  loss  Lmo  is  calculated 
using  Equation  2-15,  The  troposcatter  loss  is  calculated  in  subroutine  TROPSC. 


Table  4-1.  Summary  of  Terrain-Dependent  Propagation  Formulas  for  LOS  or  Diffraction 


Region 

Excess  Propagation  Loss 

Limits 

Equation  in 
Text 

Line-of-sight 

Alos  ”  ^ 

min  (Aseloss>  Aref) 

min  [oj  Aseloss  (Loi)  AgEF>  ArefI 

^  AreF 

fMHz  <  16 

16  :SfMHz  ^  20 
fMHz  >  20 

4-8 

Aref  ”  “20  log  ( 1  Espw/Eq  I ) 

- 

4-6 

Aseloss  from  Table  3-1 

- 

Diffracting 

-^DIF.HF  “■ 

Aseloss 

SO 

VI 

4-25 

AoiF.a  “ 

aAoip  hf  +  (1  “  Of)  Adif  VHP 

16  <  fviHz  ^  20 

4-26 

AdIF,VHF  “ 

min  (Aseloss>  Arnife) 

^  Arnife 

k  ^  3  and  Aj^p  7  dB 
Imhz  ^  20 

k  <  3  or  Aj^  >  7  dB 
fMHz  >  20 

4-24 

■^sELoss  from  Table  3-1 

- 

- 

a  =  -^4  ^ 

- 

4-9 

^KNIFE  “ 

-20  log  (|Ed.f/EJ) 

- 

4-18 

Are  **“  ^Apjj. 

- 

4-23 
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The  formulas  pertaining  to  propagation  loss  estimates  over  terrain  are  semi-empirical  and  involve 
certain  assumptions  and  approximations.  The  major  empirical  points  are  summarized  here. 

•  The  terrain-reflected  energy  is  assumed  to  emanate  from  a  region  surrounding  the  terrain 
point  for  which  the  ratio  of  clearance  to  the  radius  of  the  first  Fresnel  zone  ellipsoid  is  a 
minimum.  Actually,  the  region  formed  by  the  intersection  of  the  first  Fresnel  zone  ellipsoid 
with  the  terrain  plane  contributes  to  the  reflected  energy.  For  low  antennas,  this  region  is 
virtually  all  the  terrain  between  transmitting  and  receiving  points,  including  terrain  on  both 
sides  of  the  great  circle  plane.  Such  transverse  terrain  is  not  considered  in  TIREM. 

•  The  mean  scattering  coefficient  (Equations  4-5  and  4-16)  is  empirically  determined  from 
available  measurements,  regardless  of  the  attributes  of  the  terrain  at  the  various  measurement 
locations. 

•  When  taken  into  account,  the  surface- wave  contribution  within  LOS  is  derived  from  a  flat 
earth. 

•  In  the  diffracted  region,  each  peak  is  assumed  to  be  a  ridge  and  is  treated  as  an  ideal  knife- 
edge,  regardless  of  the  appearance  of  the  peak.  For  one  knife-edge,  the  diffraction  loss 
calculated  is  rigorous.  For  more  than  one  knife-edge,  the  method  where  the  first  knife-edge 
is  treated  as  a  source  illuminating  the  second  knife-edge  is  an  expedient  approximation. 

The  case  where  the  ridges  are  of  equal  elevation  (as  would  be  the  situation  for  plain-like 
areas,  where  each  terrain  point  acts  as  a  ridge  and  contributes  about  6  dB  loss)  is  the  worst 
situation.  A  more  rigorous  analysis  shows  that  the  loss  contribution  diminishes  as  the 
number  of  ridges  increases.'*'^ 


'••^Millington,  G.,  Hewitt,  R.,  and  Immirzi,  F.  S.,  "Double  Knife-Edge  Diffraction  in  Field-Strength 
Predictions,"  lEE  Monograph.  No.  5073,  March  1962,  pp.  419-429. 

^■^Pogorzelski,  R.  J.,  "A  Note  on  Some  Common  Diffraction  Link  Loss  Models,"  Radio  Science, 
Vol.  17,  No.  6,  Nov.  -  Dec.  1982,  pp.  1536-1540. 

'•'Wogler,  L.  E.,  "An  Attenuation  Function  for  Multiple  Knife-edge  Diffraction,"  Radio  Science, 
Vol.  17,  No.  6,  1982,  pp.  1541-1546. 


4-16 


ECAC-HDBK-93-076 


•  The  loss  at  lower  frequencies,  <  16  MHz,  is  derived  from  a  spherical  earth  regardless  of  the 
terrain  features. 

4.4  PROPAGATION  OVER  LAND/SEA 

Although  TIREM  is  intended  primarily  for  evaluation  of  the  path  loss  over  rough  earth  surfaces, 
there  may  be  situations  where  path  loss  estimates  are  needed  for  paths  that  consist  of  alternate 
segments  of  land  and  sea.  For  example,  a  path  across  a  bay  or  along  a  chain  of  islands  may  require 
the  calculation  and  combining  of  diffraction  losses  due  to  both  knife-edge  and  spherical-earth 
diffraction. 

Whenever  any  segment  of  a  propagation  path  passes  over  sea  water,  the  land/sea  propagation 
loss  is  estimated.  The  land/sea  loss  calculation  is  an  algorithm  that  provides  a  weighted  combination 
of  the  rough-earth  formulation,  described  in  paragraph  4.3,  with  the  spherical -earth  method,  taken 
from  Chapter  3.  The  presence  of  sea  water  in  the  path  profile  is  recognized  during  the  profile 
examination  process  by  observing  that  points  having  zero  elevations  are  sea  water.  Although  it  is 
realized  that  occasional  points  with  zero  elevations  may  be  present  when  a  profile  extends  into 
regions  of  below-sea-level  land,  a  sequence  of  profile  points  must  represent  at  least  1  percent  of  the 
path  before  it  is  regarded  as  sea. 

Propagation  over  land/sea  implies  that  the  propagation  path  traverses  both  irregular  terrain  over 
land  and  the  spherical-earth  surface  representing  the  sea;  antennas  may  be  on  land  or  at  sea. 
Millington'*"’®  formulated  a  semi-heuristic  method  based  on  modeling  both  land  and  sea  as  a  spherical 
earth.  Millington’s  method  successfully  predicts  the  so-called  recovery  effect  that  is  observed  when 
the  received  field  strength  increases  as  the  antenna  moves  from  land  to  sea.  This  effect  is  shown  in 
Figure  4-3  by  measurements  for  a  land/ sea  path  at  a  low  frequency  of  3.13  MHz.  Measurements 
made  over  land  are  shown  with  the  symbol  +  and  over-sea  measurements  are  shown  as  o  or  x . 

For  irregular  terrain,  Millington’s  technique  does  not  apply  and  a  practical  field  determination  is 


'‘■’“Millington,  M.  A.,  and  Isted,  G.  A.,  "Ground  Wave  Propagation  Over  An  Inhomogeneous  Smooth 
Earth,  Part  2:  Experimental  Evidence  and  Practical  Implications,"  Proc.  lEE,  Vol.  97,  Part  III, 
No.  48,  July  1950,  pp.  209-217. 
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made  in  TIREM  using  weighting  factors.  The  use  of  weighting  factors  is  described  in  paragraph 
4.4.1.  Both  techniques  predict  a  recovery  effect  and  are  compared  in  Figure  4-3  for  Millington’s 
measured  data. 


SLOUGH  (o)  NCWHAVCN 

-  - 


(6) 


OltPPl 


Figure  4-3.  Comparison  of  Millington’s  Method  and  Results  Derived  Using  Weighting  Factors 

With  Measurements  (From  Reference  4-10) 


4-18 


ECAC-HDBK-93-076 


4.4.1  Weighting  Factors 

Whenever  any  part  of  a  propagation  path  passes  over  sea  water,  the  land/sea  propagation  loss  is 
estimated  as  a  weighted  combination  of  the  propagation  loss  for  the  irregular  terrain  over  land  and 
the  spherical-earth  propagation  loss  for  the  sea.  The  propagation  loss  for  land  segments  is 
calculated  by  the  appropriate  formula  from  Table  4-1,  and  the  loss  for  sea  segments  is  taken  from 
Table  3-1.  Land/sea  paths  are  subdivided  into  three  basic  path  segments  as  will  be  discussed  in 
paragraphs  4.4.2. 1,  4. 4. 2. 2,  and  4.4. 2. 3. 

For  each  path  segment,  the  land/sea  weighting  factors  I  and  s  are  the  ratios: 


I  = 


(4-28) 


(4-29) 


so  that 


5  +  /  =  1 


(4-30) 


and  where  dLANo  and  dsEA  are  the  sums  of  the  great-circle  distances  over  land  and  sea,  respectively, 
and  d  is  the  great-circle  distance  of  the  path  segment. 

The  land/sea  weighting  factors  for  the  total  path  consisting  of  all  path  segments  are  similarly 
defined,  with  L,  S,  D,  Dland?  and  Djea  replacing  /,  s,  d,  dLANo^  and  djEA)  respectively,  in  Equations 
4-28  through  4-30. 
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Also  required  are  weighting  factors  p  and  q  for  a  transition  from  the  sea-diffracted  path  section 
to  a  LOS  path  section,  given  as: 


p  =  1  -  q 


(4-31) 


q  = 


d„ 

2d-d„ 


(4-32) 


where  do  is  the  great-circle  distance  of  the  sea-diffracted  arc  of  the  segment,  and  d  is  again  the 
great-circle  distance  of  the  path  segment.  The  coefficients,  p  and  q,  serve  to  impose  a  smooth 
transition  between  spherical-earth  LOS  and  diffraction  losses  as  the  direct  ray  just  grazes  the  surface 
of  the  sea. 

4.4.2  Propagation  Formulas  for  Land/Sea  Profiles 

The  land/sea  losses  are  heuristically  formulated  and  are  expressed  as  excess  propagation  losses 
relative  to  free  space.  The  presence  of  sea  water  in  the  terrain  profile  is  recognized  by  zero 
elevations.  Losses  for  path  segments  over  sea  water  are  calculated  using  the  spherical-earth 
formulation  discussed  in  Chapter  3.  For  BLOS  propagation  paths  leading  over  islands  with 
diffracting  peaks,  the  losses  are  calculated  in  analogy  to  the  Epstein-Peterson  method  in  that  the  total 
loss  is  the  sum  of  the  partial  diffracting  losses.  The  partial  losses  are  spherical-earth  and  knife-edge 
losses  of  a  particular  segment  of  a  propagation  path.  For  paths  undergoing  more  than  one 
diffraction,  the  spherical-earth  propagation  loss  between  the  transmitting  and  receiving  point  is  also 
calculated  and  compared  with  the  loss  derived  from  the  extended  Epstein-Peterson  method  to  select 
the  minimum  path  loss. 

A  number  of  frequently  recurring  quantities  are  required  to  describe  the  propagation 
formulations;  they  are  defined  in  Table  4-2.  The  excess  loss  calculations  for  a  given  frequency  are 
dependent  on  the  elevation,  antenna  height,  and  distance.  The  land/sea  loss  representations  are,  for 
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example,  of  the  form  Ajea  (e,  D),  implying  that  the  spherical-earth  propagation  loss  is  evaluated 
with  the  dielectric  constant  and  conductivity  of  sea  water,  with  both  antenna  heights  e  (=  elevation 
+  actual  antenna  height)  and  with  D  (=  distance)  as  parameters.  The  formulations  for  the  over-sea 
cases  are  given  below.  These  formulations  are  implemented  in  subroutine  TIRLOS  for  LOS 
land/ sea  paths  and  in  subroutine  DIF  for  BLOS  land/ sea  paths. 


Table  4-2.  Land/Sea  Definitions 


AseLOSs’ 

Spherical-earth  loss  above  free  space  with  given  e  and  a  for  land 

Asea* 

Spherical-earth  loss  above  free  space  with 
€  =  81,  (T  =  4  S/m 

A,: 

Reflection  loss  at  grazing  incidence 

(Equation  4-22) 

AF,k’ 

Knife-edge  diffraction  loss  only  (Equation  4-21) 

Ai,k’ 

Reflection  loss  (Equation  4-19  for  k  knife  edges) 

h: 

Antenna  height 

e: 

Elevation  above  MSL  -f-  antenna  height 

d,  D: 

Great  circle  distances* 

d,  D 

Linear  distances* 

/,  L: 

Land  weighting  factors  (Equation  4-28)* 

s,  S: 

Sea  weighting  factor  (Equation  4-29)* 

p,  q: 

Weighting  factors  (Equations  4-31  and  4-32) 

*Lower-case  letters  apply  to  path  segments,  capital  letters  apply  to  entire  path.  || 
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4.4.2. 1  Line-of-Sight  Land/Sea  Path 

In  the  following  figures,  the  terminating  antennas  are  represented  by  darkened  circles  and  the 
diffracting  points  are  represented  by  open  circles.  For  a  LOS  land/sea  path  (where  the  path  segment 
represents  the  total  path,  Figure  4-4),  the  excess  propagation  loss  Al/s,los  is  a  weighted  combination 
of  the  spherical-earth  loss  over  a  smooth  sea  surface  and  the  appropriate  loss  for  propagation  over 
the  land. 


■'*^L/S,LOS  ^  -^SEA  ^  -^LOS 


(4-33) 


Figure  4-4.  Line-of-Sight  Land/Sea  Path 

If  the  frequency  is  below  16  MHz,  the  spherical-earth  model  is  applicable  for  both  the  sea  and  land 
portions  of  the  path.  Therefore,  in  this  case  Al/s,los  is  simply  the  weighted  combination  of 
spherical-earth  losses  for  an  all-sea  and  an  all-land  path.  Ground  constants  appropriate  for  sea  water 
and  antenna  heights  above  sea  level  are  used  when  the  loss  is  calculated  for  an  all-sea  path. 

However,  the  user-specified  ground  constants  for  land  and  antenna  heights  measured  above  ground 
level  are  used  when  the  loss  is  calculated  for  an  all-land  path. 

■^L/s,Los  “  ^  -^sEA  ^  -^sELoss  f  <  16  MHz  (4-34) 
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On  the  other  hand,  if  the  frequency  is  above  20  MHz,  the  over-land  contribution  is  simply  the 
reflection  loss  Aref  as  calculated  by  Equation  4-6,  and  the  LOS  land/sea  loss  is: 

=  S  A,ea  +  L  (h,D)  f  >  20  MHz  (4-35) 

At  grazing  distance  D  =  Dqr,  where  by  Equation  4-7  Aref  =  A^,,  the  above  equation  becomes: 

Ae,s.los  =  S  <e.D„,)  .  L  A„  f  >  20  MHz  (4-36) 

For  frequencies  between  16  and  20  MHz  inclusive,  the  over-land  loss  is  interpolated  by 
frequency  using  the  spherical-earth  and  reflection  loss  algorithms.  The  interpolation  coefficient  o:  is 
calculated  as  in  Equation  4-9. 

■^L/s.Los  “  SA5g^(e,D)  +  L  [o:  Ajrloss^^’^^  ■^ref  (4-37) 

16  <  f  <  20  MHz 

4.4.2.2  Land/Sea  Path  Diffracted  Over  Sea  Only 

For  a  land/sea  path  with  the  diffracting  part  of  the  path  over  sea  only  (see  Figure  4-5),  the 
formulation  is  identical  with  Equation  4-34  when  the  frequency  is  below  16  MHz: 


A 


L/S.SELOSS 


SA,„  (e,D)  .  (h.D)  f  <  16  MHz 


(4-38) 


Figure  4-5.  Land/Sea  Path  Diffracted  Over  Sea  Only 
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If  the  frequency  is  above  20  MHz,  the  formulation  is: 

=  S  A^(e,D)  +  L  [p  A.  f  q  A,^(h,D)]  f  >  20  MHz  (4-39) 

At  grazing  distance  D  =  Dqr,  where  =  0  and  hence  p  =  1,  q  =  0  the  above  equation 
becomes: 

=  S  A,„(e,D^  ^  L  A„  f  >  20  MHz  (4-40) 


Equation  4-40  is  identical  with  Equation  4-36,  so  that  there  is  a  smooth  transition  between  LOS  and 
diffraction. 

For  frequencies  between  16  and  20  MHz  inclusive,  the  over-land  portion  of  the  loss  is 
interpolated  by  frequency  using  the  spherical-earth  and  reflection  loss  algorithms.  The  interpolation 
coefficient  a  is  again  calculated  as  in  Equation  4-9. 


^L/S,SELOSS  ~  ^SEA.(®»^)  ^  [(P<^  ■'■Q)^SELOSs(^’^) 

16  <  f  <  20  MHz 


(4-41) 


Equation  4-37  also  becomes  identical  to  the  above  equation  at  grazing  distance  when  p  =  1, 
q  =  0,  and  Aref  =  A,,.  Furthermore,  Equation  4-41  reverts  to  Equation  4-38  when  f  =  16  MHz 
(a  =  1)  and  to  Equation  4-39  when  f  =  20  MHz  (a  =  0).  These  methods  assure  a  smooth  loss 
transition  between  LOS  and  diffraction  over  sea  and  between  appropriate  loss  algorithms  in  the 
transition  frequency  range. 
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4.4.2.3  T.and/Sea  Path  With  Knife-Edge  and  Sea  Diffraction 

All  other  land/sea  path  configurations  are  obtained  by  cascading  the  path  segments  displayed  in 
Figure  4-6  with  land  segments  having  no  over- water  diffraction  in  various  combinations.  These 
paths  always  consist  of  a  sea-diffracted  path  segment  followed  (or  preceded)  by  at  least  one  knife- 
edge  diffraction.  The  partial  losses  of  these  subpath  segments  are  added  in  a  manner  analogous  to 
the  extended  Epstein-Peterson  method  for  knife  edges  over  land.  Let  this  excess  loss  (relative  to 
free  space)  be  designated  Amked-  The  total  excess  path  loss  is  then  expressed  as: 

■^L/S, MIXED  ~  ^  -^SELOSS^^’^^l’  -^MIXEd) 


The  term  in  the  square  brackets  represents  the  estimated  loss  of  the  total  path  over  a  spherical 
earth,  weighted  by  the  ratios  L  and  S.  This  term  tends  to  predominate  at  HE  frequencies  and 
below,  when  the  antennas  are  low  to  the  ground,  or  the  path  passes  predominately  over  sea. 

The  term  Amked  also  estimates  the  loss  for  the  total  path  but  combines  the  losses  attributable  to 
a  mixture  of  path  segments  where  either  knife-edge  diffraction  over  land  or  spherical-earth 
diffraction  over  sea  occurs.  This  mixed  loss  tends  to  predominate  when  the  surface  wave  terms  of 
the  spherical-earth  loss  become  negligible  such  as  for  high  antennas,  essentially  over-land  paths,  or 
frequencies  at  VHP  or  above.  The  BLOS  path  is  divided  into  segments,  each  of  which  includes 
either  a  single  over-sea  diffraction,  or  one  or  more  land  kmfe-edges  not  interrupted  by  over  sea- 
diffraction.  The  end  points  for  each  segment  are  defined  by  the  land  knife  edges  adjacent  to  each 
arc  of  sea  water  where  diffraction  takes  place  (see  Figure  4-6  for  examples). 

Three  cases  of  paths  involving  over-sea  diffraction  are  examined  by  TIREM.  These  are: 

•  Case  1  -  both  transmitter  and  the  receiver  horizons  are  land  (both  end  segments  are  land) 

•  Case  2  -  only  one  end-segment  of  the  total  path  is  diffracting  over  sea 

•  Case  3  -  both  end-segments  of  the  total  path  are  diffracting  over  sea. 
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The  first  case  is  simply  an  extension  of  the  method  described  in  paragraph  4.2.2  for  diffraction 
over  an  all-land  path  with  k  knife  edges.  The  method  sums  the  individual  losses  for  the  k  knife 
edges  using  the  extended  Epstein-Peterson  method  and  accounts  for  the  reflection  losses  in  the 
foregrounds  of  the  transmitting  and  receiving  antennas  by  adding  the  greater  of  these  two  losses 
(Equation  4-18).  This  procedure  is  further  extended  for  land/sea  paths  to  incorporate  any  additional 
losses  due  to  segments  with  spherical-earth  diffraction  that  may  be  distributed  among  the  knife 
edges.  Each  such  subsegment  is  bounded  by  its  adjacent  knife  edges  and  appears  as  in  Figure  4-6, 
example  A.  The  added  loss  for  each  such  subsegment  Asub  is: 

^SUB  "  ^  ^SEA  (e,D)  +  / 

where  the  weighting  coefficients  j  and  Z  are  the  respective  fractions  of  sea  and  land  along  the  profile 
of  each  such  subsegment.  Amked  for  this  first  case  is  then  simply  the  combined  knife-edge 
diffraction  losses  Akj^jp^  plus  the  sum  of  the  Agyg  losses. 

=  "^ax  (Aj  „  Aj ,)  +  5^  Ap^^  +  X;  Asub  (4-44) 

k 

If  one  or  both  end  segments  involve  diffraction  over  sea  (cases  2  or  3),  these  end  segments  are 
treated  much  like  the  land/sea  diffraction  path  discussed  in  paragraph  4.4.2.2.  This  treatment  of  the 
end  segments  having  spherical-earth  diffraction  is  necessary  to  ensure  a  smooth  transition  between 
LOS  and  diffraction,  to  employ  an  appropriate  reflection  term,  and  to  define  the  geometry  for  the 
knife  edge  that  bounds  such  a  segment. 

For  example,  if  only  one  end  segment  is  diffracting  over  sea  (case  2),  Amked  is  given  by  the 
sum  of  the  excess  losses  for  knife-edge  diffraction  (including  the  first  knife  edge  that  defines  the 
limit  of  the  end  segment),  the  sum  of  the  Ajub  losses  for  any  sea-diffracted  subsegments  along  the 
path,  and  a  term  Aqne  that  accounts  for  reflection  and  spherical-earth  diffraction  within  the  end 
segment. 


A 


MIXED 


[Ap  j  (d,dj  +  E  Ap  J  +  AguB  +  A^j^g 

k>2 


(4-45) 
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In  the  above  equation: 


Aone  =  min  {[Aj,  (d)  +  p  AJ,  Aj 

+  q[-y  Asea  (e,d)  +  /  ^SELOSS  (h,d)] 

estimates  the  losses  over  the  sea-diffracting  end  segment.  The  term  in  the  square  brackets  of 
Equation  4-45  is  the  sum  of  all  the  knife-edge  losses  for  the  path.  The  first  term,  for  the  first  knife 
edge,  shows  explicitly  the  distances  d  and  d^  that  are  illustrated  in  Figure  4-6,  example  B.  The 
distance  can  be  the  distance  to  the  next  knife  edge  or  the  distance  to  the  terminating  point.  For 
grazing  incidence,  where  p  =  1,  q  =  0,  Aqne  reduces  to  A^  so  that: 


Amked  ^Ap  k  A,^ 
k 


(4-47) 


The  result  is  the  same  as  that  obtained  when  the  over-sea  path  is  considered  as  an  all-land  path  at 
grazing.  (When  "d  of  Figure  4-6,  example  A,  is  above  grazing,  the  configuration  is  evaluated  as  an 
all-land  path  with  at  least  one  knife-edge  diffraction). 


If  both  end  segments  are  diffracting  over  sea  (case  3),  Amked  is  given  by: 


•^MDCED 


^Ap  j.  +  S  Ajue  +  A.rwo 

k 


(4-48) 
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Ap  k  is  again  the  diffracting  loss  due  to  the  k-th  knife  edge,  and  Ajwo  is: 


Atwo  =  {(Pi  P2)  K’  AJ 

+  qj  [5j  Agg^  (Spdj)  +  /j  AggLoss 
+  [ij  Agg^  (62’^  h  '^sELoss 

The  indices  1  and  2  on  p,  q,  s,  /,  e  and  d  refer  to  the  geometry  of  the  two  end  segments  each  of  the 
type  displayed  in  Figure  4-6,  example  B. 

The  example  shown  in  Figure  4-7  demonstrates  the  cascading  of  the  subpath  losses.  Two 
diffracting  peaks  (k  =  2)  are  between  two  terminating  over-sea  diffracting  paths.  There  is  no  over¬ 
sea  diffraction  between  two  diffracting  knife  edges  (i.e.,  E  Asub  =  0).  The  sum  of  the  partial  loss 
is: 

Amked  =  \,i  ^  Ag2  +  min  {(Pi  +  P2)  Aj,  Aj 

+  q^  [5j  Agg^  (Oj,  dj)  +  /j  AgggQgg  (hj,  dj)]  (4  50) 

^2  [-^2  Asea  ^2)  ^2  AgELOSS  (h^,  <i2)] 

If  there  were  an  over-sea  diffraction  between  the  two  peaks  shown  in  Figure  4-7,  then  the  Agua  term 
expressed  by  Equation  4-43  would  be  added  to  the  above  equation. 
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CHAPTER  5 

MODEL  LIMITATIONS,  STATISTICS,  MEASUREMENT 
ERROR,  AND  REPEATABILITY  STATISTICS 

5.1  MODEL  LIMITATIONS 

Propagation  models  inevitably  make  certain  assumptions  and  approximations  that  ignore  or 
simplify  the  complex  description  of  real-world  variables.  The  SEM  and  TIREM  models  are 
intended  for  estimating  path  loss  between  antennas  situated  on  a  smooth  spherical  earth  or  an 
irregular  earth,  respectively.  The  effects  of  man-made  or  natural  obstructions  other  than  the  terrain 
itself  are  not  accounted  for  (unless  the  terrain  profile  has  been  adjusted  to  account  for  these 
features),  nor  are  the  effects  of  terrain  not  located  along  the  great-circle  path  between  the  antennas. 

The  existence  of  an  atmosphere  whose  refractivity  gradient  decreases  linearly  with  height 
(Reference  1-4,  p.  4-1)  is  assumed  for  heights  below  3200  m  and  an  exponentially  decreasing 
refractivity  is  assumed  above  3200  m  (Reference  2-1,  p.  1).  For  frequencies  above  10  GHz,  the 
atmospheric  absorption  due  to  molecular  oxygen  and  water  vapor  is  included  and  is  estimated  using 
the  1962  US  Standard  atmosphere  (Reference  2-2)  adjusted  for  the  user-specified  surface  humidity. 

None  of  the  following  factors  that  may  affect  path  loss  in  practical  situations  are  considered  by 
either  TIREM  or  SEM: 

•  Ionospheric  effects 

•  Ducting  phenomena 

•  Long-term  path  loss  variability  (although  these  effects  may  be  incorporated  within  the  driver 
program) 

•  Multipath  effects  (except  in  SEM  when  the  theoretical  loss  option  for  spherical  earth  is 
selected) 

•  Absorption  due  to  rain,  foliage,  or  man-made  obstructions. 
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5.2  SEM  MODEL  STATISTICS 

Within  EC  AC,  the  standard  of  comparison  for  spherical-earth  computations  is  the  EFFSECC 
model  because  it  uses  exact  methods  for  the  loss  calculations.  Statistical  results  comparing 
propagation-loss  values  obtained  with  SEM  and  EFFSECC  for  1069  LOS  and  BLOS  paths  are 
shown  below.  This  subset  was  derived  by  stepping  the  frequency  fMHz»  the  antenna  heights  of 
transmitter  h^  and  receiver  h^ and  distance  d,,,„  in  approximately  logarithmic  steps  as  follows: 

W  =  1,  10,  10^  10^  10^ 

h,„h,,,  =  1,5,  10,  10^5X10M0' 

dk„  =  6,  12,  25,  50,  10^  2xl0^  4xl0^  8x10^,  10^ 

The  total  set  comprises  1620  paths,  from  which  paths  resulting  in  troposcatter  were  eliminated. 

For  each  of  the  1069  paths,  the  loss  difference  Ljem  -  Leffsecc  was  computed.  The  statistics  of 
the  loss  difference  are: 


PL  =  -1.3  dB,  a  =  3.5  dB 


where  n  is  the  mean  error,  and  o  is  the  standard  deviation  of  error.  The  advantage  of  SEM  is  that 
the  runtime  is  about  1/35  of  the  runtime  of  EFFSECC. 

5.3  TIREM  MODEL  STATISTICS 

TIREM  model  statistics  were  determined  by  comparing  TIREM  predictions  with  an  extensive  set 
of  propagation  data,  collected  in  various  propagation-measurement  programs.  Table  5-1  provides  a 
summary  of  the  different  propagation-measurement  programs  including  the  location,  the  terrain 
characteristics,  the  frequencies,  path  length,  etc.  The  data  selected  for  comparison  included  data 
from  each  program  (all  available  frequencies  and  distances  are  used,  but  not  all  available  antenna 
heights  and  only  data  obtained  using  copolarized  antennas). 
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Table  5-1.  Summary  of  Terrain  Characteristics 


Location 

Reference® 

Description  of  Terrain 

Worldwide 

5-1 

No  foliage  effects  are  included  since  the  data  comprise  measurements  obtained  by 
monitoring  reliable  communications  links  over  long  (yearly)  periods. 

Path  Length:  10  -  800  km 

Frequency:  40  -  4000  MHz 

Flagstaff, 

Arizona 

5-2 

Area  is  characterized  by  extremely  rough  terrain  interspaced  with  flat  areas.  Foliage 
varies  from  dense  forest  to  widely  spaced  trees  to  sparsely  covered  land  in  flat  areas. 
Transmitter  and  receiver  sites  are  relocated. 

Path  Length:  1  -  30  km 

Frequency:  2,  4,  8,  16,  32,  64,  128,  256,  400,  512,  800  MHz 

Colorado 
(Gunbarrel  Hill) 

5-3 

Receiver  site  is  near  the  summit  of  a  hill  in  the  open  plains  of  the  Rocky  Mountain 
foothills.  Ten  of  the  55  transmitter  sites  associated  with  this  receiver  site  are  located 
in  the  mountains  and  only  one  of  these  results  in  a  line-of-sight  path.  Five 
transmitter  sites  are  concealed  by  foliage. 

Path  Length:  1  -  120  km 

Frequency:  230,  400,  800,  1800,  4600,  9200  MHz 

Colorado 
(North  Table 

Mt.) 

5-4 

Receiver  site  is  on  a  high  mesa  at  the  juncture  between  the  mountains  and  plains  and 
visible  to  most  of  its  associated  59  transmitting  sites. 

Path  Length:  1  -  120  km 

Frequency:  230,  400,  800,  1800,  4500,  9200  MHz 

Colorado 

(Plains) 

5-5 

Transmitter  site  is  on  a  plain.  Area  is  composed  of  smoothly  rolling  hills  east  of 
transmitter.  No  vegetation  or  trees. 

Receiver  sites  at  fixed  semicircular  distances,  approximately  120  sites. 

Path  Length:  5  -  80  km 

Frequency:  20,  50,  100  MHz 

Colorado 

(Mountains) 

5-6 

Same  transmitter  site  as  in  plains.  Area  is  very  rugged  mountains  west  of 
transmitter.  Woods  and  trees  present,  approximately  40  sites. 

Path  Length:  5  -  50  km 

Frequency:  20,  50,  100  MHz 

Ohio 

5-7 

Irregular  and  partly  wooded  terrain.  Part  of  measurements  taken  with  snow  on 
ground.  Receiver  sites  in  approximate  semicircles,  with  five  peripheral  transmitters. 
Approximately  250  sites. 

Path  Length:  5  -  50  km 

Frequency:  20,  50,  100  MHz 

^References  5-1  through  5-7  are  provided  at  the  conclusion  of  Table  5-1. 
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Table  5-1.  Summary  of  Terrain  Characteristics  -  Continued 
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The  results  of  the  comparison  of  TIREM-predicted  losses  with  measurements  are  summarized  in 
Table  5-2  for  each  of  the  propagation  modes.  The  statistics  shown  in  this  table  are  the  mean  and 
the  standard  deviation  of  the  loss  difference,  Ltirem  -  LMeasurements-  For  all  the  modes,  the  overall  loss 
difference  had  a  mean  of  -0.6  and  standard  deviation  of  10.5  dB. 

In  Table  5-2,  mean  and  standard  deviation  for  the  spherical-earth  propagation  mode  and  for 
f^Hz  ^  16  are  displayed  separately  to  demonstrate  that,  for  lower  frequencies,  the  spherical-earth 
approximation  is  good  (all  the  measurements  are  in  the  Flagstaff  area  characterized  by  rough 
terrain). 


Table  5-2.  TIREM  Statistical  Summary 


Propagation 

Mode 

Nominal 

Frequency 

Number 

of 

Measurements 

Loss  Difference 

Mean 

dB 

Standard 

Deviation 

dB 

Spherical  Earth 

2,  4,  8,  16 

102 

-2.2 

6.9 

Line-of-sight 

20,  32,  50,  64,  100, 

230  400,  800,  1800 
4600,  9200 

1217 

-2.8 

8.9 

Diffraction 

2798 

0.2 

11.4 

Troposcatter 

60  -  4000 

358 

1.5 

8.8 

Total 

2  -  9200 

4475 

-0.6 

10.5 

5.4  ERROR  ESTIMATES  OF  THE  MEASUREMENTS 

The  Bell  Aerosy stem’s  report  on  propagation-loss  measurements  includes  a  detailed  analysis  of 
the  instrumentation  errors  associated  with  transmitting  and  measuring  systems.  The  reported 
accuracies  of  path-loss  measurements  as  associated  with  instrumentation  errors  were  shown  to  be 
within  ±4.1  dB  for  frequencies  from  2  MHz  to  32  MHz  and  within  ±3.2  dB  for  frequencies  from 
32  MHz  to  900  MHz.  This  report  clearly  conveys  the  difficulties  inherent  in  the  error  analysis  of 
large-scale  propagation-measurement  programs,  and  is  the  only  measurement  analysis  readily 
available. 


5-5 


ECAC-HDBK-93-076 


5.5  STATISTICS  CONCERNING  REPEATED  MEASUREMENTS 

The  Bell  Aerosystem’s  report  also  contains,  for  each  path,  a  series  of  repeated  measurements. 
Measurement  equipment  was  returned  to  the  previous  sites.  For  southern  Arizona,  repeated 
seasonal  measurements  were  taken  at  approximately  3-month  intervals.  In  northern  Arizona,  a 
series  of  measurements  was  taken  every  5  weeks.  Measurements  of  the  seasonal  variations  were 
taken  less  frequently  in  southern  Arizona  because  this  region  has  just  two  temperature  seasons  per 
year,  which  vary  only  slightly;  whereas,  in  northern  Arizona  the  four  seasons  are  more  varied,  with 
cool  summer  months  and  considerable  snow  in  winter. 

The  results  of  the  repeated  measurement  series  for  frequencies  from  2  to  800  MHz  are  shown  in 
Table  5-3.  It  shows  that  the  repeatability  is  worse  when  the  frequency  increases. 
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APPENDIX  A 

INPUT  PARAMETER  RANGES 

A.l  INPUT  PARAMETER  RANGES 


The  allowable  ranges  for  the  input  parameters  for  TIREM  are  shown  in  Table  A-1,  and  for  SEM 
in  Table  A-2. 


Table  A-1.  TIREM  Input  Parameter  Ranges 


Variable 

CONDUC 


EXTNSN 


HPRFL 


Conductivity  of  the  earth 
surface 


Description 


Valid  Range 
0.00001  to  100  S/m 


Profile  indicator  flag: 

.TRUE.  -  profile  is  an  extension  of  the  last  path  along  a  radial 


.TRUE,  or  .FALSE. 


.FALSE.-  new  profile 

Array  of  profile  terrain 
heights  above  mean  sea  level 


-450  to  9000  m 


HUMID 

Surface  humidity  at  the  transmitter  site 

0  to  50  g/m^ 

NPRFL 

Total  number  of  profile  points  for  the  entire  path 

>3 

PERMIT 

Relative  permittivity 

1  to  100 

POLARZ 

Transmitter  antenna  polarization: 

'V'  or  'H' 

PROPFQ 

RANTHT 

REFRAC 

TANTHT 

XPRFL 


"V"  -  vertical 

'H"  -  Horizontal  _ 

Transmitter  frequency  _ 

Receiver  structural  antenna  height _ 

Surface  refractivity _ _ 

Transmitter  structural  antenna  height _ 

Array  of  great-circle  distances  from  the  beginning  of  the  profile 
point  to  each  profile  point.  _ _ 


1  to  20000  MHz 
>0  to  30000  m 
200  to  450  N-units 
>0  to  30000  m 
>  0  m 
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Table  A-2.  SEM  Input  Parameter  Range 


Variable 

Description 

Valid  Range 

CONDUC 

Conductivity  of  the  earth  surface 

0.00001  to  100  S/m 

HUMID 

Surface  humidity  at  the  transmitter  site 

0  to  50  g/m^ 

LANDHT 

Average  terrain  elevation  of  the  spherical  earth  above  mean  sea 
level 

0  to  9000  m 

PERMIT 

Relative  permittivity 

1  to  100 

POLARZ 

Transmitter  antenna  polarization: 

'V'  -  vertical 
'W  -  horizontal 

'V'  or 

PROPFQ 

Transmitter  frequeney 

1  to  20000  MHz 

PTHLEN 

Great-circle  distance  from  transmitter  to  receiver 

>  0  m 

RANTHT 

Receiver  structural  antenna  height 

>0  to  30000  m 

REFRAC 

Surface  refractivity 

200  to  450  N-units 

SEOPTN 

Spherical-earth  output  option: 

1  -  free-space  wave  loss 

2  -  envelope  of  the  field  maxima  (loss  minima) 

3  -  envelope  of  the  field  minima  (loss  maxima) 

4  -  theoretical  loss  for  spherical  earth 

1  to  4 

j 

TANTHT 

Transmitter  structural  antenna  height 

>0  to  30000  m 

A-2 


ECAC-HDBK-93-076 


APPENDIX  B 

TROPOSPHERIC  SCATTER 


B.l  GENERAL 

The  long-term,  median,  basic  transmission  loss  for  paths  where  forward  scatter  in  the 
troposphere  is  the  dominant  coupling  mechanism  is  given  by  Equation  2-15. 


=  max  {[30  log  f  -  20  log  D  +  F(D0)  +  20  log  (rJD)],  LpJ  +  H„  -  F„  (B-1) 


The  empirical  functions  F(D0),  Hq,  and  Fo  are  described  in  this  appendix. 

B.2  ATTENUATION  FUNCTION  F(D0) 

The  attenuation  function  F(D0)  is  an  empirical  function  of  the  product  of  the  mean  sea  level  arc 
distance  D  and  the  scatter  angle  Figures  B-1,  B-2,  B-3,  and  B-4  contain  plots  of  the  attenuation 
function  F(D0)  and  the  coefficients  for  the  values  of  =  250,  301,  350,  and  400.  The 
empirical  curves  are  approximated  by  a  fourth-order  polynomial  as: 


¥(06)  =  [\n(DdJ10)T 

k=o 


(B-2) 


Evaluation  of  ¥(06)  for  intermediate  values  of  and  s  is  by  linear  interpolation;  for  s  >  1  the 
value  of  s  is  replaced  by  1/s.  The  path  asymmetry  parameter  is: 


s 


oc. 


(B-3) 
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s 

bo 

bl 

b2 

b3 

0.01 

172.1598 

15. 178922 

0.673849 

-0.02665 

0.0065 

0.02 

172. 0214 

15.095866 

0.748526 

0.01245 

0.0045 

0.05 

172.0847 

15.610355 

0.788888 

-0. 162934 

0.0721 

0.10 

172.3286 

16.509260 

0.752966 

-0.42048 

0.1741 

0.20 

172.65585 

16.661147 

0.719552 

-0.310936 

0.2228 

0.30 

172.62052 

16.116216 

0.851304 

*0.04552 

0.1817 

0.50 

172.99412 

16.805433 

0.756862 

-0.089631 

0.2832 

0.70 

173.08063 

17.334882 

0.796767 

-0.246248 

0.3435 

1.00 

173.05149 

17.192346 

0.823508 

-0.162147 

0.3401 

Figure  B-1.  Attenuation  Function,  F  (D0),  and  the  Coefficients  for 
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s 

t>0 

b2 

“’3 

t>U 

0.01 

168.57 

15.569 

0.6M67 

-0.17929 

0.014284668 

0.02 

168,61 

15.285 

0.5828 

0.029297 

0.0157242 

0.05 

168.37 

15.557 

1 .01042 

-0.069279 

0.04924011 

0.10 

168.10 

16.525 

1.62945 

-0.528870 

0.1628723 

0.20 

168.39 

17.253 

1.7478 

-0.6342316 

0.2445526 

0.30 

168.70 

17.938 

1.8719 

-0.693817 

0.2944489 

0.50 

168.78 

18.130 

1.8984 

-0.671310 

0.3582764 

0.70 

168.85 

18.5H3 

1.9567 

-0.837936 

0.4245453 

1.00 

168.91 

18.666 

1.9751 

-0.855865 

0.4409485 

Figure  B-2.  Attenuation  Function,  F  (D6),  and  the  Coefficients  \  for  =  301 


s 


16*1.27 

16*1.37 

16*1.68 

163.9*1 

164.95 

164.91 

164.77 
164.73 

164.78 


15.387 

15.139 

15.247 

16.283 

16.286 

17.406 

19.547 

19.994 

20.256 


0.7816 

0.5234 

0.4200 

2.3161 

1.1963 

2.1280 

3.5410 

3.8334 

4.0000 


-0.21782 

0.08553 

0.3346 

-0.7292 

0.42007 

-0.23667 


-2.0837 


0.0556259 

0.017525 

-0.000305 

0.18669 

0.0368042 

0.182373 

0.500091 

0.619080 

0.648529 


Figure  B-3.  Attenuat 
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where 

j  =  O';  +  AQ!;  (B-4) 

The  correction  coefficient  Aa;  is  required  to  account  for  the  effects  of  a  nonlinear  refractivity 
gradient.  Aa^  is  approximated  for  different  ranges  of  j  as  follows: 

Range:  0<^„;<0.1 


Ao-j  =  max  {d^,  j(B  +  Cd^,  j)  -  A,  0} 


(B-5) 


Where  A,  B,  and  C  are  empirical  functions  of  i  as  plotted  in  Figure  B-5  for  =  250,  301, 
350,  400.  (These  plots  are  converted  into  tabular  format,  and  intermediate  values  are  retrieved 
by  linear  interpolation). 

Range:  d^,<0 

The  same  format  as  Equation  B-5  is  used,  replacing: 


^C,i  ^C.i  I  ^e^o.i 


(B-6) 


and  with  A,  B,  and  C  evaluated  at  do  -,  =  0 
Range:  0.1  <  6^-,  <0.9 

Aaj  =  =  0.1)  (B-7) 

+  N^(9.97  -  cot  0^  i)[l  -  exp(0.05  d^  j]  x  lO’® 
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Figure  B-5.  Plots  of  the  Functions  A,  B,  and  C 
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B.3  SCATTERING  EFFICIENCY  CORRECTION  F„ 

The  scattering  efficiency  correction,  Fq,  which  allows  for  the  reduction  of  scattering  efficiency  at 
great  heights  in  the  atmosphere,  is  given  by: 

F.  =  1.086  (VhJ  (h„  -  h,  -  h„  -  h,,)  (B-8) 

where 

=  0.5696  hjl  +  (0.031  -  2.32  X  lO'^  (B-9) 

+  5.67  N/  X  10-^)e-3*  h/  x  lO'®] 

and 

(B-10) 

h„  -  s  D0/(1  +  s)^ 

hjj  =  sdj0/(l  +  s)^  (B-11) 

ds  =  D  -  -  d^,2  (B-12) 

B.4  FREQUENCY-GAIN  FUNCTION  H„ 

The  frequency-gain  function  is  a  function  of  the  path  asymmetry  parameter  s  (defined  by 
Equation  B-3),  (defmed  by  Equation  B-9)  and  the  parameters  V;  and  q,  where: 

Vj  =  47r  hj  aj\  (B-13) 

-  V,/V, 


q 


(B-14) 
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with  hi  the  structural  antenna  height  of  the  transmitters  and  receivers,  respectively,  and  X  is  the 
wavelength. 

Different  methods  are  used  to  calculate  depending  on  jj^,  s,  and  q.  The  various  methods  are 
given  below. 

For  TT.  >  1  and  so  ^1.  Ho  is  given  by: 

=  [H„(Vj)  +  H„(V,)]/2  +  AH„  (B-15) 

with 

AH„  =  6(0.6  -  log  Tj)  (log  s)  (log  q),  (B-16) 


The  empirical  frequency-gain  function  Ho(Vi)  is  plotted  in  Figure  B-6  for  =  1,  2,  4,  5,  10, 
15,  20,  30,  50,  100  and  each  plot  is  approximated  by  a  fourth-order  polynomial  for  intermediate 
values  of  Vj  and  by  its  asymptotes  for  very  large  or  very  small  values  of  Vj. 


Ho(Vi)  = 


0, 


E  d,  [In(Vi)]^ 


V.  >  v„ 


,  V,  <  V,  <  V, 


k=0 


-40  log  V;  +  aj, 


V.  <  V, 


(B-17) 
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The  coefficients  dk  are  tabulated  as  part  of  Figure  B-6.  The  coefficients  aj  and  the  upper  and 
lower  limits  of  Vj  are  tabulated  in  Table  B-1.  For  intermediate  values  of  17,  Ho(Vi)  is  linearly 
interpolated  between  the  curves  of  Figure  B-6. 


Table  B-1.  Parameters  for  Calculating  Ho  For  Very  Large  and  Very  Small  Values  of  V; 

Using  Equation  B-17  


Vs 

ai 

V, 

Vu 

1 

-  0.1705530256 

0.018 

9.0 

2 

5.8108634830 

0.0255 

10.0 

4 

11.9285966381 

0.037 

12.0 

5 

13.6631059879 

0.042 

13.0 

10 

21.9423214234 

0.066 

18.0 

15 

27.6657765715 

0.09 

25.0 

20 

30.3110521994 

0.115 

40.0 

30 

36.2956340498 

0.16 

70.0 

50 

44.8655247182 

0.25 

90.0 

100 

56.9580753795 

0.48 

100.0 

The  following  constraints  apply  to  Equation  B-15: 


1. 

s 

=  min  (s,  10) 

2. 

q 

=  min  (q,  10) 

3. 

s 

=  max  (s,  0.1) 

4. 

q 

=  max  (q,  0. 1) 

5. 

If 

AHo  >  [Ho  (Vi)  -H 

Ho 

(V2)]/2  and  AHo  >  0, 

then  Ho  =  H^  (Vj)  + 

Ho 

(V2) 

6. 

If 

lAHol  <  [Ho(Vi) 

+ 

Ho  (V2)]/2  and  AHo  <  0, 

then  Hq  =  0 
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For  Ti.  >  1  and  so  =  1 .  the  functions  H(,(Vi)  are  approximated  by  a  fifth  order  polynomial  for 
intermediate  values  of  Vj  and  by  its  asymptotes  for  very  large  or  small  values  of  Vj. 


H„(Vi)  = 


0, 

5 


V.  >  Y, 


£  d,  [ln(V,)]‘  ,  V,  <  V.  <  V, 


k=0 


-40  log  Vj  -  a^, 


V.  <  V, 


(B-18) 


The  plots  of  the  function  Ho(Vi)  and  the  tabulation  of  the  coefficient  d^  for  ri^  -  1,  2,  4,  6,  10,  15 
each  for  s  =  0.10,  0.25,  0.50,  0.75,  1.00  are  found  in  Figures  B-7  through  B-12. 

For  values  Vj  that  are  very  large  or  very  small,  Ho(Vi)  is  evaluated  by  using  the  parameters 
tabulated  in  Tables  B-2  and  B-3  for  the  above  and  s  values.  For  intermediate  values  of  rj^  and  s, 
Ho(Vi)  is  derived  by  a  linear  interpolation  for  each  variable. 

For  n,  <  1  and  SQ  1 .  Ho  is  calculated  by: 

H„(,.  <  1)  =  .  0)  +  vM(t,  ■  I)  -  =  0)]  (B-19) 


where  Ho(i7s  =  1)  is  determined  from  Equation  B-15  and  where: 


H„(r,3  =0)  =  10  log  {2(1  -  s2q2)/[r2'(h(rj)  - 


(B-20) 


with 


rj  =  V,(l  +  1/s) 


(B-21) 


r^  =  VjCl  +  s)  =  qsr, 


(B-22) 
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s 

do 

d2 

d3 

d4 

‘‘5 

0.10 

11.25 

-18.99859 

9.16108 

-0.05206 

-0.77583 

0.08075 

0.25 

12.65 

-23.28189 

13.07615 

-0.09284 

-2.13960 

0.50549 

0.50 

15.35 

-27.78663 

12.92337 

2.43753 

-2.89846 

0.46697 

0.75 

17.75 

-30.07131 

11.50693 

3.36909 

-2.11519 

0.08582 

1.0 

19.ii5 

-32.12501 

11.97889 

4.05866 

-2.77689 

0.25943 

Figure  B-10.  Frequency-Gain  Function,  H,,,  for  sq 


1,  175  =  6 
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-zb.eissi 

-30.81(175 

-35.15910 

-36.997K 

-37.K0957 


7.0740(( 

10.42451 

7.79812 

6.85575 

3.43747 


2.76938 

3.90438 

5.56981 

6.15626 

6.76033 


-0.33396 

-2.47659 

-0.38235 

-0.08357 

1.63919 


<^5 

-0.35943 

0.22601 

-0.75589 

-0.89875 

-1.53728 

4 

hI! 


iS 


b9! 

K8I 


in 

H- .  i 


1  ■ 
L  1 

■  (iH 

!■ 

i 

Figure  B-12.  Frequency-Gain  Function,  for  sq 
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Table  B-2.  Upper  Limit  of  V;  for  Calculating  Ho  for  Very  Large  Values  of  Vj 

(Equation  B-17  for  sq  =  1) 


s  =  0.1 

s  =  0.25 

s  =  0.50 

s  =  0.75 

s  =  1.0 

V, 

Vu 

Vu 

Vu 

Vu 

Vu 

1 

1.9 

3.8 

7.0 

12.0 

15.0 

2 

5.0 

7.0 

10.0 

15.0 

17.0 

4 

15.0 

15.0 

15.0 

15.0 

18.0 

6 

18.0 

18.0 

18.0 

18.0 

20.0 

10 

24.0 

24.0 

24.0 

24.0 

35.0 

Table  B-3.  Lower  Limit  of  V;  and  Parameter  &2  for  Calculating  Hq  for  Very  Large 


Values  of  V;  (Equation  B-17  for  sq  =  1) 


Vs 

■ 

■ 

■ 

0.25 

a2 

■ 

■ 

■ 

■ 

1 

0.1 

16.25 

0.1 

10.75 

0.1 

5.45 

0.122 

1.78 

0.14 

1.25 

2 

0.1 

10.75 

0.1 

5.25 

0.14 

7.20 

0.172 

4.29 

0.2 

7,15 

4 

0.1 

4.75 

0.149 

1.83 

0.209 

7.35 

0.25 

10.86 

0.295 

18.80 

6 

0.13 

6.52 

0.20 

6.79 

0.265 

11.79 

0.32 

15.12 

0.37 

17.86 

10 

0.18 

5.27 

0.27 

12.08 

0.38 

17.95 

0.46 

21.44 

0.53 

24.03 
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The  function  h(ri)  is  given  in  the  form  of  a  plot  as  shown  in  Figure  B-13  and  is  approximated  by: 


where 


h(r.)  = 


1.630637  r^,  ri<  0.001 

exp(-0.47  +  0.445X  -  0.1152x^ 

-  0.007954x^),  0.001  <ri< 0.6887 

exp{  -exp(v)/1000} ,  0.6887  ^  r^  <  100 


exp  ( - 1 . 906295/r  1 ) , 


r;>100 


X  =  In(rj) 


(B-23) 

(B-24) 


=  6.18705  -  0.892717X  -  0.284649x^  +  0.023584x^ 


(B-25) 


For  1 .  SQ  =  1 .  and  s  ^  1 .  Hg  is  calculated  as  in  Equation  B-19,  except  Hg(Tj5  0)  is  given 
by: 


=  0)  =  11.68  -  7.792X  +  1.310x^  +  0.04222x3H„(7/^)  =  -0.016545x'‘  (B-26) 

where 

X  =  ln(r,) 


For  r).  <  1.  sq  =  1  and  s  >  1.  Hg  is  calculated  as  above  in  Equation  B-26,  with: 

X  =  InCr^) 


(B-28) 
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APPENDIX  C 

SEM  CRITICAL  DISTANCES 


C.l  GENERAL 

Within  LOS,  the  space  wave,  as  given  by  Equation  3-1,  oscillates  with  distance  (see  Figure  3-2). 
In  SEM,  there  is  an  option  (SEOPTN  =  4)  to  calculate  the  path  loss  corresponding  to  the  theoretical 
space-wave  lobing  structure  Espw/E,,  that  results  from  alternate  constructive  and  destructive 
interference  between  the  direct  and  reflected  waves  within  the  LOS  region  over  a  perfectly  smooth 
spherical  earth. 

When  this  option  is  chosen.  Equation  3-1  is  used  to  calculate  Egpw/Eo  for  any  distance  less 
than  dspw  Alternative  choices  calculate  losses  corresponding  to  the  envelope  of  the  space-wave  field 
maxima  Emax/Eo  (SEOPTN  =  2)  or  the  envelope  of  the  space-wave  field  minima 
Emin/Eo  (SEOPTN  =  3). 

With  SEOPTN  =  2,  Equation  3-14  is  used  out  to  the  distance  dwAx  where  the  last  field 
maximum  occurs  and  Equation  3-1  is  used  between  d^Ax  dgpw 

With  SEOPTN  =  3,  Equation  3-15  is  used  out  to  the  distance  d^m  where  the  last  field  minimum 
occurs  and  Equation  3-1  is  used  between  dj^iN  and  djpw 

A  fourth  option  (SEOPTN  =  1)  permits  use  of  a  free-space  approximation  Epsw/Eo  to  the 
oscillating  space-wave  field  within  the  reflection  region  to  a  maximum  distance,  doNE,  where  the 
space-wave  field  last  becomes  equal  to  the  free-space  field.  Equation  3-1  is  again  used  beyond  doNE 
to  the  distance  dspw 

This  appendix  outlines  the  procedures  for  finding  the  maximum  range  at  which  the  normalized 
space- wave  electric  field  has  a  value  of  1  (D  =  doNE).  the  range  of  the  last  maximum  (D  =  dwAx)  of 
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the  normalized  field,  and  the  range  of  its  last  minimum  (D  =  dwiN)-  The  distance  d^^x  is  also  used 
in  the  calculation  of  djp^,  given  by  Equation  3-17. 

C.2  QUARTIC  EQUATION  METHOD 

This  approximate  method  for  finding  d^Ax  is  sufficiently  accurate  for  situations  involving 
propagation  over  land  and  for  frequencies  generally  above  HE  or  for  horizontal  polarization.  It  is 
used  wherever  possible  because  the  calculations  are  much  faster  than  the  more-accurate  iterative 
method,  which  is  required  for  propagation  of  vertically-polarized  waves  over  sea  at  lower 
frequencies. 

For  antennas  near  the  ground  and  for  distances  that  are  great  compared  to  the  antenna  heights, 
such  that  r/  ~  rj',  the  phase  factor  A  (Equation  3-4)  may  be  approximated  by: 


where 


A  = 


4  7rh/h/ 
X  D 


d 


MAX 


d,  +d 


2 


d,/d2  =  hj/h^ 


(C-1) 


(C-2) 


(C-3) 


(C-4) 


(C-5) 
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From  Equation  3-1  it  is  evident  that,  for  plane  earth  and  assuming  a  reflection  coefficient  of 
R  =  -1,  the  last  oscillation  maximum  occurs  at  A  =  tt,  so  that: 


4  X  h/  h/ 


X  d, 


=  TT 


MAX 


(C-6) 


Using  Equations  C-2  through  C-5,  Equation  C-6  can  be  brought  into  the  form: 


^max  “  2  C  d^iAx 


Xa,(hi  +h,) 


Cd„*x  =0  (C-7) 


where 


^  ^  (hj  +  h^)^  (C-8) 

hA 

The  Quartic  Equation  C-7  is  then  solved  using  Brown’s  method  of  biquadratic  equations. The 
solution  of  a  quartic  equation  yields  four  (4)  roots,  of  which  all  four  (4)  roots  are  complex,  or  two 
(2)  roots  are  complex  and  two  (2)  roots  are  real,  or  all  four  (4)  roots  are  real.  Since  the  solution  is 
a  physical  quantity  (i.e.  distance),  it  is  assumed  that  real  roots  will  always  exist.  Also  the  solution 
must  be  greater  then  zero  but  less  then  the  radio  horizon  distance,  therefore,  the  distance  to  the  last 
oscillation  peak  is  the  largest  real  root  less  then  the  radio  horizon  distance.  This  has  proved 
empirically  correct. 


<^-'Brown,  W.  S.,  "Solution  of  Biquadratic  Equations,"  Aircraft  Engineering,  January,  1944,  page  14. 
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The  last  oscillation  minimum  occurs  near  A  =  Itt,  which  yields  the  following  quartic  equation: 


•*MAX 


-  2Cd 


MAX 


2Xa^(h,  +112) 


-4a,(h,.h2)C  =0 


(C-9) 


Since  the  distance  to  the  last  oscillation  minimum  must  be  greater  then  zero  and  less  then  the 
distance  to  the  last  oscillation  maximum,  the  largest  real  root  satisfying  this  condition  is  used. 

C.3  ITERATIVE  METHOD 

For  some  cases  of  propagation  of  vertically-polarized  waves,  the  assumption  that  the  specular 
reflection  coefficient  is  -1  can  be  invalid.  For  these  cases,  the  quartic  solution  can  result  in  an 
erroneous  estimate  for  d^Ax-  Incorrect  values  of  d^Ax*  itt  him,  can  cause  the  path  loss  predicted  by 
SEM  to  exhibit  inconsistent  behavior.  For  these  cases,  an  efficient,  iterative  algorithm  is  used  to 
calculate  the  range  parameters  of  interest.  The  algorithm  requires  more  time  resources  than  the 
quartic  solution  method,  but  yields  the  correct  answer  in  all  cases. 

To  minimize  the  number  of  times  that  the  iterative  method  is  used,  the  quartic  solution  method  is 
first  called  to  make  an  initial  range  estimate  for  all  cases. 

For  vertical  polarization,  this  range  estimate  is  checked  for  consistency  with  the  expected 
characteristics  of  the  normalized  space  wave  field  at  its  last  maximum  (the  magnitude  of  the  last 
maximum  >  1  for  example)  and  the  iterative  algorithm  is  called  if  a  contradiction  is  detected. 

First,  a  method  to  find  the  greatest  range  at  which  the  normalized  electric  field  has  a  value  equal 
to  1,  doNE>  is  described.  Then  this  value  is  used  to  find  the  range  to  the  last  maximum,  d^Ax^ 
the  range  to  the  last  minimum,  dwiN- 


C-4 


ECAC-HDBK-93-076 


Last  Range  At  Which  |  E/En  |  =  1 

This  section  describes  an  iterative  procedure  for  finding  the  greatest  range,  doNE*  at  which  the 
normalized  magnitude  of  the  electric  field  is  equal  to  1.  Using  the  normalized  field  values  shown  in 
Figure  3-2  as  typical,  it  is  recognized  that  doNE  occurs  at  some  distance  less  than  the  radio  horizon 
given  by  Equation  3-18. 

The  following  information  is  used  to  locate  the  relevant  distance  parameters. 

•  The  field  variation  in  the  region  between  doNE  and  dnoR  is  nearly  linear  on  a  linear/log 
scale 

•  The  radian  path  length  difference,  27r(r2’  -  r/l/X,  monotonically  decreases  to  zero  at  the 
horizon 

•  The  phase  of  the  reflection  coefficient  changes  slowly  as  a  function  of  distance  in  the 
vicinity  of  d^Ax 

•  If  the  change  of  phase  of  the  reflection  coefficient  between  dnoR  and  d^iN  is  small,  the 
radian  path  length  difference  between  the  radio  horizon  and  d^iN  is  about  lir 

The  distance,  doNE.  is  found  by  iteratively  decreasing  the  distance  from  the  horizon  until  the 
normalized  space-wave  field  |  Espw/Eq  |  is  1.  Because  |  Esp^/Eq  |  oscillates  through  1  many 
times,  its  derivative  with  respect  to  distance,  D,  and  the  phase  difference  A  between  the  estimated 
distance  of  doNE  and  dnoR  are  kept  within  specified  limits  to  ensure  that  doNE  is  the  distance  nearest 
the  horizon  where  |  Espw/Eq  |  last  becomes  1 . 

The  normalized  space-wave  field  is  calculated  from  Equation  3-1,  but  the  effects  of  the  path 
length  ratio  and  divergence  factor  are  ignored  (the  product  (ri7r2’)Dc  =  1)  since  they  affect  only  the 
magnitude  of  the  field  and  not  the  range  locations  of  interest. 
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To  avoid  problems  with  numerical  precision  in  the  estimation  of  the  value  of  the  derivative  (with 
respect  to  distance)  of  the  normalized  field  at  dnoR,  a  starting  distance  of: 


■•start 


ad 


HOR 


(C-10) 


is  used,  where  a  is  initially  0.975.  It  is  important  that  dgxART  be  between  doNE  and  dnoR,  and  far 
enough  from  dnoR  that  the  iteration  will  converge.  To  this  end,  the  value  of  a  is  adjusted  until  the 
phase  A  is  between  ir/20  and  ir/2  and  the  normalized  field  is  between  0.1  and  1. 

Next,  a  range  step  dr  is  determined  according  to: 


dr 


=  d 


START 


|fl^|)/slope  *  log 


(C-11) 


where 

(C-12) 

The  estimated  maximum  distance  djEsi  where  the  normalized  space-wave  field  last  becomes  1 
within  the  horizon  is  then: 

d  =  d  -  dr  (C"13) 

'^lEST  ^START 


slope  = 


1 


log  (1  1/dg.j.^j) 


dD  'E„ 
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Equations  C-11  and  C-12  are  derived  assuming  a  linear  variation  in  the  normalized  field  with  the  log 
of  distance  between  dnoR  and  doNE-  Because  the  field  variation  is  not  exactly  linear,  Equation  C-13 
will  sometimes  yield  a  distance  less  than  the  true  value  of  doNE-  This  is  because  the  variation  of  the 
field  is  not  exactly  as  assumed. 

To  avoid  overshooting  the  true  value  of  doNE>  a  range  step  limit  is  calculated  based  upon  the 
change  of  A  between  dsxART  and  djEsr- 

Assuming  that  the  path  length  difference  between  the  direct  and  reflected  fields  is  the  dominant 
factor  influencing  the  change  of  phase  of  |  E/Eq  |  near  doNE.  then  the  resultant  field  will  increase  to 
a  value  of  1  as  distance  decreases  below  dnoR  accompanied  by  a  change  in  phase  of  less  than  tt. 

The  total  change  of  phase  is  defined  as: 


Ay  A(dj£5y)  A(djjQjj) 


(C-14) 


where  A  =  27r(r2’  -  ri’)/X  as  in  Equation  3-4.  The  phase  difference  Ay  is  not  permitted  to  exceed 
TT.  To  this  end,  if  |  Ay  [  >  tt  the  range  step  is  limited  as  follows: 


—  -^(^HOR  tiiEsr) 

Ay 


(C-15) 


The  range  step  is  limited  to  the  minimum  of  Equations  C-11  and  C-15.  Of  course,  if  the  step  is  so 
limited  then  Ay  must  be  recalculated  and  checked  that  its  value  is  less  than  tt.  Although  the  above 
limit  on  the  range  step  may  delay  the  achievement  of  convergence,  it  ensures  that  the  distance  doNE 
will  not  be  bypassed. 
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The  iteration  proceeds  by  substituting  the  last,  best  estimate  of  doNE  for  ^start  iii  Equation  C-11 
to  find  a  range  step,  limit  the  step  if  necessary,  and  calculate  the  value  of  the  normalized  space- 
wave  field  as  described  above.  The  iteration  stops  when: 


|.(d  =  -  1|  <  0.01  (C-16) 

^0 


The  above  procedure  is  implemented  in  the  FORTRAN  routine  RICALC.  This  routine  contains 
additional  logic  that  is  necessary  to  handle  certain  unique  space-wave  field  distributions  (limiting 
situations  such  as  antenna  heights  equal  to  zero,  frequency  limits,  etc.). 

Range  of  Last  Maximum 

The  last  maximum  of  the  normalized  space-wave  field  |  Egpw/Eo  |  occurs  at  the  distance  d^Ax 
when  the  direct  and  reflected  fields  are  in  phase  for  the  last  time  before  the  antennas  become 
separated  beyond  the  horizon.  At  d^Ax*  the  value  of  |  Esp^/Eg  |  lies  between  1  and  2  (the 
maximum  it  could  be  if  the  earth  were  a  perfect  reflector).  Also,  at  dMAx>  h  is  known  that  the  first 
derivative  of  the  field  must  change  signs.  Therefore,  the  search  for  d^Ax  begins  at  doNE  by  stepping 
back  in  range  (with  an  appropriate  range  step  limit)  until  a  change  in  the  sign  of  the  first  derivative 
of  the  field  is  located.  A  bisection  algorithm  is  used  to  refine  the  calculation  of  d^Ax-  this  case, 
the  starting  point  is: 


d 


START 


=  d, 


ONE 


(C-17) 


and  the  range  step  dr  is: 


dr 


=  d 


START 


I  |)/slope  +  log 

Eo 


(C-18) 
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and  d^iEST* 


^MEST  ^ONE 


(C-19) 


is  an  estimate  of  d^Ax- 

A  criterion  is  needed  to  determine  if  the  range  step  is  reasonable.  Again,  this  criterion  is 
determined  from  available  information  about  the  change  in  phase  of  the  field  over  the  range  step: 


^(‘^MESt)  ^(^ONe)  ^  ^ 

where  Aj  is  the  total  phase  change  between  doNE  and  the  estimate  for  the  range  of  the  last  maximum 
dMEST-  Examination  of  a  typical  space-wave  field  distribution,  such  as  in  Figure  3-2,  shows  that  the 
difference  between  the  phase  of  the  field  at  doNE  and  dwAx  is  less  than  or  equal  to  tt. 

Therefore,  a  range  step  limit: 


^  ^(^one  ‘^mest) 

At 


(C-21) 


is  imposed. 

The  approximation  dwEsi  is  a  reasonable  estimate  of  the  range  of  the  last  maximum  dwAx-  To 
check  the  estimate,  the  value  of  the  derivative  of  the  field  at  the  ranges  d^Ax  and  doNE  are  calculated 
and  compared  for  a  sign  difference.  If  there  is  a  sign  difference,  then  the  location  of  the  maximum 
has  been  successfully  bracketed  and  a  bisection  algorithm  can  be  used  to  resolve  the  location  of  d^Ax 
to  any  desired  accuracy.  If  there  is  no  sign  difference  between  the  values  of  the  derivatives  at  the 
two  range  points,  then  it  is  assumed  that  another  range  step  is  necessary.  Another  range  step  is  then 
calculated  using  Equations  C-19,  C-20,  and  C-21  with  the  current  estimate  of  d^Esr  substituted  for 
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doNE  to  obtain  a  new  estimate,  and  the  above  procedure  is  repeated  until  a  range  bracket  for  d^Ax  ts 
found. 

The  calculation  of  the  range  of  the  last  field  maximum  is  implemented  in  the  FORTRAN  routine 
REXCAL.  This  routine  also  contains  additional  logic  to  handle  exceptional  cases. 

Range  of  Last  Minimum 

Using  the  same  procedure  that  led  to  finding  the  range  of  the  last  field  maximum,  we  may  find 
the  range  of  the  last  field  minimum.  To  ensure  that  the  search  begins  at  a  range  point  on  the 
transmitter  side  of  the  last  field  maximum  (which  is  necessary  to  use  the  same  algorithm),  begin  the 
search  at: 

d  =  0  9d  (C-22) 

’^START 


The  minirmim  normalized  field  will  never  be  less  than  0  so  a  range  step  is  calculated  by: 


dr 


=  d 


START 


D/slope  +  log  (d„^x) 


(C-23) 


The  estimate  of  the  distance  to  the  last  field  minimum  dwEST  is: 


^mEST  ^START 


dr 


(C-24) 


and  the  range  step  is  limited  by: 


dr  < 


IT 


(d: 


START  '^mEST^ 


(C-25) 
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Using  this  new  estimate  for  the  range  of  the  last  field  minimum,  the  derivative  of  the  field  is 
calculated.  The  derivative  at  the  initial  range  (dsTARx)  is  also  calculated  and  the  signs  of  the  two 
derivatives  are  compared.  If  a  sign  change  is  detected,  a  bisection  algorithm  is  used  to  refine  the 
estimate.  If  no  sign  change  is  observed  between  the  two  values  of  the  derivative,  the  above 
procedure  is  repeated  until  the  required  sign  change  indicating  the  field  minimum  is  observed. 

The  calculation  of  the  range  of  the  last  field  minimum  is  implemented  in  the  FORTRAN  routine 
REXCAL.  An  argument  is  passed  through  the  parameter  list  to  indicate  whether  the  range  of  the 
last  maximum  or  the  range  of  the  last  minimum  is  desired. 
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APPENDIX  D 

REFLECTION  LOSS 


D.l  PHASE  FACTOR 


The  phase  difference  b  between  the  direct  wave  and  the  wave  reflected  from  a  plane  surface  is 
(Figure  D-1): 


6  =  ^  [(r^  +  r,)  -  (dj  +  d^)]  (radians) 


(D-1) 


where 


d,  +  d^  =  d 


(D-2) 


= 


d,  .  d.  >  >  c 


(D-3) 


r 


2 


>  >  C 


(D-4) 


so  that 


(D-5) 
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Figure  D-1.  Reflection  by  Plane  Surface 


The  width  f  of  the  first  Fresnel  zone  is  defined  as  the  ground  clearance  c  for  which  the  phase 
factor  is  6  =  tt.  Therefore,  from  Equation  D-5. 


(D-6) 


and  the  phase  difference  can  be  written  in  terms  of  the  ratio  of  the  clearance  c  of  the  direct  ray 
above  the  reflecting  point  to  the  first  Fresnel  zone  width  f. 


d  =  7r(c/f)^ 


(D-7) 
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For  di  >  >  hi  and  d2  >  >  h2,  the  distances  dy  and  d2  can  be  replaced  by  Dj  and  D2  (see  Figure 
D-1),  so  that  Di  +  D2  =  D.  The  clearance  c  and  the  width  of  the  first  Fresnel  zone  f  can  then  be 

expressed  as; 


c 


hi  + 


h2  hj 


D, 


(D-8) 


f 


X/ 

1 

--1 

N 

Di 

D2 

(D-9) 


The  phase  difference  5  has  a  minimum  at  6  =  A 


A  =  TT  min(c/f)^ 


(D-10) 


for  which 


dA 

dD^ 


0 


This  condition  results  in  the  requirement  that; 


b.  =  h. 

D.  D, 


(D-11) 


(D-12) 
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This  is  the  geometrical  optical  condition  that  reflection  takes  place  when  the  angle  of  the  incident 
ray  equals  the  angle  of  the  reflected  ray.  Equation  D-10  is  therefore  the  phase  difference  for 
interference  of  the  direct  and  reflected  ray. 

D.2  Plane-Earth  Theory  Space-Wave  Reflection  Loss 

The  effect  of  ground  reflection  is  approximated  in  TIREM  by  assuming  that  a  single  dominant 
reflecting  plane  exists  at  a  terrain  profile  point  that  has  the  smallest  ratio  of  direct-path  clearance  to 
first  Fresnel  zone  width,  c/f.  As  shown  above  in  D.l,  the  phase  difference  is: 

A  =  TT  min(c/f)^  (D-13) 

The  resulting  reflection  loss  Arep  is  calculated  using  Equation  4-6  in  which  the  normalized 
space- wave  field  |Espw/Eo|  can  be  calculated  using  either  Equation  4-2a  or  Equation  4-2b.  When 
Equation  4-2a  is  used,  the  field  oscillations  that  occur,  as  the  clearance  ratio  c/f  changes,  become 
evident  because  of  the  cosA  term  in  that  equation.  The  reflection  factor  -Arep  is  plotted  from 
Equation  4-2a  in  Figure  D-2  using  dashed  lines  for  Po  =  ^  and  Pg  =  0.2636.  These  values  of  Pg  are 
the  limiting  values  for  frequencies  of  1  MHz  and  20  GHz,  respectively,  calculated  with  Equation 
4-5.  Note  that  the  Pg  =  1  curve  is  identical  with  the  plane-earth  theory  reflection  factor  curve  for 
R  =  1  published  by  Bullington.'^  '  ®'^ 

Because  the  field  oscillations  that  occur  in  real  irregular  terrain  are  not  predictable  with  accuracy 
using  plane-earth  reflection  theory,  they  should  be  eliminated.  Elimination  of  these  oscillations  can 
be  accomplished  by  substituting  the  free-space  loss  for  cases  where  the  clearance  ratio  exceeds 
VTTT.  This  value  of  c/f  yields  a  phase  difference  A  of  7r/3  which,  when  substituted  in  Equation 
4-2a  with  Pg  =  I,  results  in  a  normalized  space- wave  field  of  1  (i.e.,  a  free-space  field).  This 


'^■'Bullington,  K.,  "Radio  Propagation  Fundamentals,"  The  Bell  System  Technical  Journal.  Vol.  36, 
No.  3,  May  1957,  p.  608. 

^^Bullington,  K. ,  "Radio  Propagation  for  Vehicular  Communications, "  IEEE  Transactions  on  Vehicular 
Technology.  Vol.  VT-26,  No.  4,  November  1977,  p.  303. 
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substitution  is  indicated  in  Figure  D-2  by  the  dark  line  extending  to  the  right  of  c/f  0.577 
(Fresnel  zone  1/3)  at  the  ordinate  -Arep  =  0. 

In  TIRFM,  the  free-space  field  is  substituted  for  clearance  ratios  c/f  ^  VT7¥  =  0.564.  This 
value  is  used  because  the  approximate  Equation  4-2b  is  actually  used  in  TIRFM  to  estimate  the 
reflection  loss  from  free-space  for  low  clearance  ratios.  Substitution  ofVTTlr  in  Equation  D-13 
results  in  A  =  1  and,  with  Po  =  1,  Equation  4-2b  then  predicts  that  the  normalized  space  wave  field 
is  equal  to  the  free-space  field.  Figure  D-2  shows  curves  of  -Aref>  calculated  using  Equation  4-2b, 
for  selected  frequencies  from  1  MHz  to  20  GHz.  These  are  the  solid  lines  extending  to  the  left  of 
the  point  c/f  =  0.564,  Aref  =  0.  It  is  seen  that  TIREM  includes  Aref  only  when  the  clearance 
ratio  is  less  than  vTTlr  and  that  the  scattering  coefficient  p^  diminishes  the  effect  of  ground 
reflection  as  the  frequency  increases. 
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